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I. 


INTRODUCTION 


For  the  past  year  Mission  Research  Corporation  has  been  pursuing 
a  computational  and  theoretical  investigation  of  high  power  microwave 
interactions  with  air  for  the  Naval  Research  Laboratory.  The  technical 
objectives  were  to  design,  construct,  operate  and  benchmark  against 
experimental  data  a  one-dimensional,  self-consistent  computer  simulation 
of  microwave  breakdown  of  nitrogen.  A  self-consistent  nitrogen 
chemistry-wave  optics  hydrodynamic  Microwave  Gas  Breakdown  simulation  code 
MINI  has  been  developed.  We  have  used  MINI  to  simulate  the  MRC 
experiments  being  performed  for  NRL  (ref.  1)  and  also  to  study  the 
microwave-plasma  coupling  from  atmospheric  to  low  (25  Torr)  pressure. 

Preliminary  simulations  have  been  performed  with  a  simplified 
version  of  MINI  which  excludes  hydrodynamic  motion  of  the  electrons  and 
the  neutral  gas.  The  qualitative  features  of  these  "No-Hydro"  simulations 
appear  to  be  consistent  with  experimental  results.  These  preliminary 
simulations  are  extremely  useful,  since  they  allow  both  hydrodynamic 
effects  and  effects  due  to  plasma  reflectivity  to  be  easily  isolated.  The 
No-Hydro  simulations  were  used  primarily  to  study  the  reflection  and 
absorption  of  the  microwaves  with  the  breakdown  gas.  This  is  an 
approximation  even  for  short  pulses  since  the  neutral  gas  requires  only  a 
small  drift  speed  to  conduct /convect  large  amounts  of  energy  out  of  a 
region.  In  principle  MINI  could  be  run  for  very  long  times  to  follow  the 
neutral  gas  motion;  however,  times  much  larger  than  1  usee  become  costly. 

In  the  course  of  investigation,  the  electron  density  profile  was 
seen  to  be  reasonably  modeled  using  either  a  symmetric  or  slab  model 
profile.  Using  a  slab  model  for  the  electron  density  profile  allowed  a 
simple  simulation  to  be  performed  using  a  separately  developed  long-term 
chemistry  code.  The  issue  addressed  in  these  long-term  chemistry 
simulations  (LTCS)  was  the  reflectivity  of  the  breakdown  plasma  due  to 
thermal  ionization  of  the  neutral  gas.  In  a  thermodynamic  equilibrium, 
with  a  gas  temperature  of  one  electron  volt,  the  breakdown  plasma  would 


reflect  essentially  all  of  the  incident  microwave  energy.  These 
simulations  were  typically  run  for  tens  of  microseconds  and  provided 
valuable  quantitative  Insight  into  gas  heating  by  microwave  radiation.  In 
what  follows,  a  brief  synopsis  of  the  code  development  and  the  theory  to 
date  on  microwave  gas  breakdown  will  be  discussed. 

II.  COOE  DEVELOPMENT 

The  primary  objective  in  fiscal  year  1981  was  the  development  of 
a  self-consistent  nitrogen  chemistry  -  wave  optics  -  hydrodynamic 
simulation  to  investigate  microwave  decoupling  due  to  gas  motion. 

However,  research  was  initiated  using  a  wave  optics  chemistry  code  to 
address  the  issue  of  plasma  reflectance  from  a  high  density  plasma.  The 
results  showed  that  the  energy  deposition  of  the  microwave  could  not  be 
sufficiently  localized  to  heat  the  gas.  The  details  of  the  wave  optics 
chemistry  code  are  given  in  Appendix  A.  The  code  follows  the  electric 
field  by  solving  the  wave  equation  for  a  plane  wave  of  normal  incidence. 
The  electromagnetic  boundary  conditions  were  of  the  Lindman  (Ref.  2) 
type.  The  electric  field  provides  the  energy  source  for  the  11  species 
nitrogen  chemistry  code,  CHEMN2,  developed  for  this  project.  CHEW12 
calculates  the  time  development  of  ground  state  molecular  nitrogen,  atomic 
nitrogen,  their  first  charge  state,  the  excited  states,  N2+(B),  N4+, 

N(2D),  NjjCA),  NgCB),  N^C);  and  the  electrons.  The  heavy  particle 
temperature,  Tg,  vibrational  temperature,  Tv,  and  the  electron 

temperature,  Te>  are  a^so  calculated.  The  set  of  stiff  differential 
equations  are  solved  using  a  modified  version  of  CHEMEQ  (Ref.  3).  The 
wave  optics  code  has  been  thoroughly  tested  and  is  used  extensively  in  the 
simulation  of  the  short  pulse  experiments. 

The  development  of  the  hydrodynamic,  nitrogen  chemistry  and  wave 
optics  simulation  MINI  (Microwave  Nitrogen  Interaction  Code)  was  completed 
in  early  1982.  As  an  extension  of  the  wave  optics  chemistry  simulation 
MINI  follows  the  hydrodynamics  of  the  previous  mentioned  species. 

Transport  processes  for  the  gas  include  viscous  heating  and  thermal 
conduction.  The  transport  coefficients  were  taken  from  Bird  (Ref.  4  and 


5).  A  Euler Ian  differential  equation  of  the  hydrodynamic  equation  in 
conservative  form  is  solved  using  the  algorithm  of  Rubin  and  Burstein 
(Ref.  6).  The  scheme  is  simple,  second  order  in  accuracy,  and  is 
relatively  stable  while  centering  the  dissipation  terms.  Although  no  use 
Is  made  of  any  artificial  viscosity,  the  code  can  be  switched  to  use  the 
Fiux-Corrected-Transport  (FCT)  algorithm  of  Boris  (Ref.  7).  The 
hydrodynamic  and  chemistry  portions  are  calculated  separately  in  a  cycle 
assuming  sufficiently  weak  flows.  The  hydrodynamics  is  determined  from 
the  updated  values  of  the  chemistry.  Hence,  the  hydrodynamics  is 
considered  as  a  small  correction  to  the  chemistry. 

MINI  follows  both  the  electrons  and  gas  hydrodynamics  In  a 
two  fluid  model  assuming  separate  fluids  for  the  electrons  and  the  heavy 
particles.  The  continuity  equation  for  each  species  is  solved.  A  single 
mass-average  velocity  for  the  heavy  particles  is  calculated.  The 
diffusion  of  the  neutral  species  is  neglected  except  for  the  dominant 
molecular  nitrogen  species.  The  contribution  from  the  electric  field  is 
Included  in  describing  the  diffusion  of  the  charged  species.  The  DC 
electric  field  Is  assumed  to  be  ambipolar  with  the  sum  of  the  positively 
charged  particle  flux  equal  to  the  electron  particle  flux.  In  the 
strongly  colllsional  region,  the  electrons  do  not  sample  regions  of 
different  electric  field  or  pressure.  We  assume  the  electrons  are  in 
steady  state  as  are  the  positively  charged  species.  The  ambipolar  field 
assumption  allows  for  the  solution  of  the  ambipolar  field  and  the  electron 
drift  speed.  Sources  which  cool  the  electrons  and,  therefore,  heat  the 
positively  charged  particles  and,  hence,  the  gas  include  the  following 
terms:  Joule  heating  by  AC  and  DC  electric  fields,  elastic  and  inelastic 
heating,  and  work  due  to  momentum  transfer.  Since  the  translational  and 
rotational  degrees  of  freedom  of  the  heavy  particles  are  closely  coupled 
and  are  in  near  equilibrium,  only  an  energy  equation  describing  the 
vibrational  energy  and  translational  energy  of  the  gas  is  included. 

Hence,  the  three  energy  equations  solved  are  for  the  electrons,  heavy 
particles,  and  the  total  vibrational  energy. 
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The  AC  electric  field  can  be  solved  in  the  regular  planar 
geometry  or  In  spherical  geometry.  However,  the  hydrodvnamlcs  and 
chemistry  are  always  solved  in  a  one  dimensional  planar  geometry.  The 
solution  of  the  spherical  wave  equation  is  used  only  In  modeling  spherical 
or  focused  systems  similar  to  the  MRC/NRL  experiments.  The  solution  to 
the  hydrodynanlc  equations  in  planar  geometry  and  the  solution  to  the  wave 
equation  in  spherical  geometry  are  not  consistent.  One  can  easily  imagine 
systems  in  which  the  lateral  scale  lengths  of  the  system  are  long  compared 
to  the  longitudinal  (parallel  to  the  propagation  vector)  scale  lengths. 

In  such  systems,  diffusion  would  actually  occur  most  rapidly  in  the 
longitudinal  direction  or  radial  direction  if  calculating  the  electric 
field  in  spherical  geometry. 

HINI  can  also  be  run  using  a  dilated  electron  mass.  This  allows 
for  longer  time  steps  in  the  chemistry  and  hydrodynamics.  Dilating  the 
electron  mass  can  be  justified  physically  when  the  electron  drift  speed  is 
slow,  i.e.,  the  kinetic  energy  is  small  compared  to  the  electron 
temperature.  As  mentioned  previously,  the  hydrodynamics  is  taken  as  a 
correction  to  the  chemistry.  Hence,  in  the  chemistry  portion  of  the 
calculation,  the  real  electron  mass  is  explicitly  used. 

III.  QUALITATIVE  FEATURES  OF  GAS  BREAKDOWN 

Based  on  simulations  using  HINI,  we  have  documented  some 
qualitative  features  of  microwave  breakdown  of  gases.  The  first  of  these 
Is  microwave  energy  coupiing  to  the  breakdown  plasma.  One  general 
characteristic  of  microwave  coupling  to  a  breakdown  plasma  is  the 
inability  to  localize  the  energy  deposition.  The  problem  of  interest 
requires  the  gas  to  be  significantly  heated  before  the  microwave  energy 
decouples.  Further  gas  heating,  as  in  a  slow  combustion  wave  (in  which 
the  gas  diffusion  time  is  comparable  to  the  gas  heating  time),  is  to  be 
distinguished  from  a  fast  detonation  wave  (in  which  the  gas  diffusion 
time  is  long  compared  to  the  gas  heating  time.)  The  latter  may  have 
substantial  military  applications.  The  basic  question  addressed  was  under 


what  circumstances  can  a  large  amount  of  microwave  energy  be  deposited 
rapidly  Into  the  gas.  Under  high  Irradlance,  the  electron  density  rises 
sufficiently  rapidly  to  partially  reflect  the  incident  radiation.  The 
total  electric  field  then  proceeds  to  break  down  the  gas  at  a  position 
away  from  the  Instantaneous  peak  electron  density  In  the  direction  of  the 
{-.ource  (anti-parallel  to  the  wave  vector).  This  secondary  Ionization 
rises  rapidly  to  absorb  the  microwaves  and  shields  (decouples)  the  primary 
Ionization  f.ont  from  the  microwave  radiation.  This  then  limits  the 
deposition  time  of  the  radiation  into  the  gas.  Secondary  Ionization 
fronts  due  to  this  process  were  observed  in  preformed  plasmas  or  plasmas 
formed  in  a  highly  spherical  (focused)  system. 

The  primary  gas  heating  mechanism  at  relatively  high  pressures 
(greater  than  25  Torr)  in  nitrogen  Is  quenching  of  the  electronic  states 
by  molecular  and  atomic  nitrogen.  Hence,  rapid  gas  heating  requires  high 
electron  density  and  temperature,  sufficiently  high  to  excite  the 
electronic  states  through  electron  impact.  Full  hydrodynamic  simulations 
of  short  pulse  nitrogen  breakdown  experiments  with  no  target  have  shown 
that  the  gas  Is  not  sufficiently  heated  to  excite  these  states. 

Simulations  show  electron  temperatures  on  the  order  of  1-2  eV  and  an 
electron  density  on  the  order  of  1  x  10u  cm-'*  for  the  35  GHz  experimental 
parameters.  The  electric  field  Is  rapidly  attenuated  in  a  distance  of 
approximately  one  wavelength.  The  results  are  expanded  in  more  detail  in 
Appendix  B.  These  simulation  results  are  qualitatively  consistent  with 
the  experiment. 

Simulation  of  microwaves  Incident  on  a  reflecting  surface  shew 
breakdown  occurlng  at  the  peak  nodes  of  the  electric  field  (the 
characteristic  standing  wave  pattern),  as  would  be  expected.  As  the 
electron  density  rises,  these  "pancakes"  of  electron  density  begin  to 
attenuate  the  radiation  such  that  the  pancake  closest  to  the  reflecting 
surface  (one-quarter  wavelength  away)  is  shielded  from  the  radiation.  As 
the  absorption  continues,  the  remaining  pancakes  would  then  sequentially 
decouple  the  next  pancake  closest  to  the  surface,  and  the  process  repeats 

itself.  Again,  in  these  low  pressure  simulations,  the  gas  is  not 
significantly  heated. 
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IV. 


COMPARISON  CP  EXPERIMENT  TO  MINI 


*  Comparison  of  MINI  to  the  NRL/MRC  experiments  has  been 
undertaken  to  verify  that  MINI  includes  the  proper  physics  to  sltmlate 
microwave  breakdown  of  nitrogen.  Comparison  of  the  simulation  predictions 
and  experimental  measurements  of  electron  density  are  shown  in  Figure  1. 
This  data  is  based  on  both  3  GHz  and  35  GHz  experiments.  Agreement  Is 
good,  degrading  towards  the  high  pressure  region.  One  explanation  for 
this  discrepancy  is  that  MINI  assumes  a  Maxwellian  distribution  for  the 
electrons.  This  assumption,  particularly  for  the  tall  electrons  which  do 
the  ionization,  may  not  be  valid  for  v/w  >  1  (high  pressure). 

Experimental  observations  of  the  spectra  have  been  unable  to 
detect  the  Na*  391hA  band.  Only  the  second  positive  neutral  nitrogen 
bands  have  been  observed.  This  is  consistent  with  MINI'S  prediction  of 
I391H/I337l  *  10~H  for  the  experimental  parameters.  This  effectively  sets 
an  upper  bound  on  the  electron  temperature  in  the  experiment,  Te  _<  2  eV. 

Finally,  ionization  fronts  such  as  those  predicted  by  MINI  were 
observed.  This  is  discussed  in  detail  in  Appendix  C  where  a  copy  of  our 
upcoming  article  is  reproduced. 
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Microwave  Energy  Deposition,  Breakdown  and  Heating  of  Nitrogen 

C.  L.  Yee 

Mission  Research  Corporation 
A.  W.  All 

Naval  Reseach  Laboratory 
ABSTRACT 

A  comprehensive  code  is  developed  for  the  microwave  energy  deposition  in 
N2  to  describe  the  breakdown  and  heating  of  nitrogen.  The  microwave 

absorption  and  reflection  from  a  nitrogen  plasma  is  based  on  wave  optics. 
The  ionization  and  the  deionization  of  N2,  however,  is  described  in  detail 
by  an  N2  chemistry  code  (CHEM-N2)  which  provides  plasma  diagnostics 

capabilities  by  calculating  emission  intensities  of  the  second  positive 
and  the  first  negative  band  systems  of  N2  and  N2+. 

The  code  is  applied  to  breakdown  calculations  in  N2  and  the  results  are  in 
reasonable  accord  with  experimental  data.  The  wave  optics  code  for  energy 
deposition  is  also  coupled  to  an  air  chemistry  code  (CHMAIR)  to  describe 
the  breakdown  off  a  reflecting  surface. 

An  insight  is  provided  into  the  plasma  absorption  by  using  a  simple  slab 
density  profile  for  the  electrons.  This  models  roughly  the  profile 

expected  in  air  breakdown  off  a  reflecting  surface.  Simulation  results 

show  a  steep  localized  electron  density  profile  at  the  interference 
maximum.  These  simple  calculations  show  the  plasma  reflectivity  to  be 
small  (<  25%)  for  1  cm  light  at  high  electron  densities,  Ne  1015  cm-3  > 
Nc  =  1013  cm-3 .  Simulation  results  also  show  tne  possibility  of 
attaining  high  electron  density.  However,  in  breakdown  off  a  reflecting 
surface,  the  microwave  is  decoupled  from  the  plasma  when  the  absorption 
length  is  less  than  one-half  the  wavelength  of  the  light.  Important 
questions  to  be  addressed  in  the  future  include  the  possibility  of  plasma 
maintenance,  gas  heating,  shock  formation,  plasma  decoupling  and  effects 
of  strong  non-linear  forces. 
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MICROWAVE  ENERGY  DEPOSITION,  BREAKDOWN 
AND  HEATING  OF  NITROGEN  AND  AIR 

I.  INTRODUCTION 

The  interactions  o£  the  high  power  microwave  radiation  with  air, 
matter  and  plasmas  have  potential  for  various  applications,  and  hence  merit 
detailed  studies  and  analysis.  A  detailed  study,  theoretical  and 
experimental,  is  necessary  in  order  to  delineate  possible  regimes  for 
applications  purposes. 

Studies  can  be  performed  analytically  to  provide  answers  in  the 
limiting  cases.  However,  the  phenomena  of  the  microwave  interaction  with 
air,  matter  and  plasmas  are  complex  and  require  time  and  space  dependent 
treatments  of  a  large  number  of  physical  processes.  This  can  be  carried 
out  in  a  self  consistent  manner  by  developing  detailed  computer  codes  which 
treat  part  or  all  the  phenomena. 

This  report  deals  with  and  describes  computer  codes,  developed  at  NRL, 
for  the  analysis  of  the  pulsed  microwave  interaction  with  nitrogen  (N2) . 

The  interaction  is  considered  from  the  first  principles  with  the  microwave 
absorption  by  one  or  few  free  electrons,  present  in  the  gas,  the  electron 
build  up,  the  breakdown,  and  the  microwave  absorption  and  reflection  by  the 
plasma  electrons.  The  computer  code  which  describes  this  phenomena  is 
basically,  consists  of  two  codes  coupled  together.  One,  describes  the 
electromagnetic  radiation  (E  and  M  Code)  and  the  other  describes  the  physics 
and  chemistry  of  the  ionization  and  the  deionization  in  N2  (CHEM-N2)  Code. 
The  E  and  M  Code,  however,  can  also  be  coupled  into  an  air  chemistry  code  to 
describe  the  breakdown  in  air,  as  discussed  in  the  report. 


Manuscript  submitted  July  24, 1981. 


A-l 


In  Section  1  the  E  and  M  Code  and  the  solution  of  the  wave  equation  is 
described.  Section  2  gives  the  description  of  the  pure  N2  chemistry  code, 
and  Section  3  deals  with  the  breakdown  calculations  in  Ns,  where  the  break¬ 
down  power  thresholds  are  compared  with  experimental  data  for  several  nitro¬ 
gen  pressures.  ,  In  Section  4,  the  E  and  M  Code  coupled  to  an  air  chemistry 
code  is  utilized  to  calculate  breakdown  in  air  near  a  surface.  For  air 
breakdown  near  a  surface  one  can  utilize  microwave  powers  below  the  threshold 
for  air  breakdown  by  a  factor  of  k  and  still  be  able  to  break  down  the  air. 
This  is  due  to  the  fact  that  the  microwave  reflection  from  the  surface 
interferes  constructively  with  the  outgoing  wave  resulting  in  doubling  the 
electric  field  at  a  distance  of  from  the  surface. 

2.  The  E  and  M  CODE 

In  order  to  describe  the  microwave  absorption  by  a  plasma  one  must  know 
the  electric  field  at  all  positions  in  the  plasma.  To  do  this  we  consider 
the  solution  of  one  dimensional  wave  equation1 


+  v 


+  to* 


u)2  +  v- 


t»2  +  ve 


0  (1) 


for  a  plane  wave  incident  normal  to  a  plasma.  The  electron  neutral  (electron- 
ion)  collision  frequency  is  denoted  by  v  and  the  plasma  frequency  is  defined 

C 

as  m2*»  — - £-  •  where  the  notation  is  standard.  Since  the  absorption 
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length  at  microwave  frequency  is  long  at  the  critical  density  =  m  u)a/^e2 
we  can  specify  an  essentially  free-space  boundary  condition2  for  the 


incoming  light  as 


=  2 


where  Uin  is  the  amplitude  of  the  known  incident  light.  This  boundary 
condition  is  also  applicable  to  a  plasma-vacuum  interface  as  constructed 
in  experiments  using  a  mylar  divider.  The  right  side  boundary  condition  in 
the  plasmas  specified  as 
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wher*  J  «  cE2  «  ~5r7”  Ez*  A8ain  ^or  l°w  densities,  the  right 

hand  side  of  equation  (2)  is  similar  to  a  free-space  propagating  to  the  right. 
Equation  (1)  is  differenced  implicitly  to  retain  stability  for  large  times 
steps  in  a  high  density  plasma,  ^  a  vc2  +  uj2.  The  implicit  difference  form 
of  equation  (1)  is 
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where  (5  s  Ej+^  ”^Ej  +  Ej  1*  T^lis  can  also  be  cast  into  a 

tri-diagonal  system  of  the  form 
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and  v  *  ~g“  “fe  +  va  A  solution  of  a  tri-diagonal 


system  of  equations3 
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can  be  shown  to  be 
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where  Che  G 
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s  end  F 
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s  satisfy  the  recursion  relation 
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The  solution  technique  is  then  to  specify  the  boundary  condition  end  sveep 
the  mesh  space  for  the  G^'s  end  F^'s.  After  which  we  sweep  again  to  obtain 
the  updated  electric  field. 

The  stability  analysis  of  the  implicit  scheme  can  be  made  simple  by 
dropping  the  last  term  of  Equation  (1).  This  term  is  typically  not 
important  in  determining  the  stability  criterion.  We  rewrite  rquation  (1) 
without  this  term  a&  a  system  of  two  coupled  equations. 
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Choosing  an  arbitrary  fourier 
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where  V  = - and  W  ** 
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component  we  let  v”  *  §nexp  (-ikjAy)  and  Wj  ■  p  exp  (-ikjAy)> 
We  find  after  some  algebra  that 
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where  v  m  cAtsln  (kAy/2)  .  Stability  requires  the  eigenvalues  of  the 

Ay 

amplification  matrix  A  to  be  less  than  or  equal  to  1.  One  can  show  that  the 
eigenvalues  are  bounded  by 
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and  hence  the  implicit  scheme  is  unconditionally  stable  for  arbitrary  or. 

The  boundary  conditions  are  explicitly  differenced  and  are  coupled  to 
an  explicit  form  of  Equation  (1)  to  obtain  boundary  value.  This  explicit 
form  of  Equation  (1)  is: 
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Once  the  boundary  values  are  obtained,  we  use  the  implicit  difference 
equation  for  all  interior  points.  The  boundary  values  are  used  as  "first" 
values  in  sweeping  the  tri-diagonal  system  of  equation.  Since  locally  the 
explicit  differencing  is  stable  for 
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the  boundary  must  be  at  low  density  to  maintain  stability, 
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In  practice,  tha  Implicit  scheme  must  satisfy  the  explicit  stability  criterion 
to  maintain  accuracy  and  minimum  phase  error.  In  considering  microwave  break¬ 
down  off  a  reflecting  surface ,  this  is  important  since  the  interference 
maximum  determines  the  breakdown  field. 

As  a  check  on  the  code,  we  perform  simulations,  on  free-space,  propaga¬ 
tion  and  reflection  off  metal  surface  in  a  vacuum.  Figure  (la)  shows  the 
reflected  light  for  a  propagating  electromagnetic  wave  in  vacuum  for  a  single 
pass  to  the  reflecting  surface  and  back  to  the  point  of  entry  of  the  radia¬ 
tion.  The  system  length  is  kL  »  1$  (where  k  *  2rr/A)  *°  the  return  time  of  the 
radiation  is  tRau  *  30  for  cAt  “  0 . 9Ay .  The  code  tracks  the  propagation 
velocity  well,  however,  in  the  transient  time  u>t  £  30,  the  numerics  has  intro¬ 
duced  errors  in  the  reflected  wave.  The  average  reflected  wave  should  be 
sero  during  the  transient  period,  yet  we  see  the  average  reflected  wave  of 
roughly  10$  of  the  incident  light.  Figure  (lb)  shows  the  spitial  electric 
field  (in  normalized  variables).  Since  the  launched  incoming  wave  is  unity, 
the  combined  wave  once  reflected  from  the  surface  should  be  of  amplitude  two. 
We  see  the  implicit  scheme  has  introduced  large  errors  in  the  phase  of  the 
wave  and  therefore  the  reconstructed  wave.  Figure  (2a)  shows  a  similar  test 
with  a  reduced  time  step  cAt  «  0.2Ay.  However,  in  this  test,  the  reflected 
wave  in  the  transient  period  is  less  than  1$  of  the  incident  light.  The 
smaller  phase  error  reconstructs  the  interference  patterns  to  its  proper 
maximum  of  two  as  can  be  seen  in  Figure  (2b).  An  extensive  parameter  study 
shows  good  results  are  obtained  using  the  implicit  scheme  for  cAt  £  0.2Ay. 

Also  as  a  check,  we  consider  light  normally  incident  on  a  linear  density 

profile.  The  density  profile  in  this  collisionless  plasma  (vc  *  O.ltu)  rises 

from  zero  to  the  cut  off  density  N  with  a  density  scale  length  of  kL  »  10. 

c  s 

The  Airy  function  solution  requires  standing  wave  pattern  toward  the  Incident 
light  and  evanescent  wave  in  the  overdense  plasma.  Figure  (3)  shows  the 
solution  obtained  for  the  electric  field  using  the  implicit  scheme  with  the 
boundary  conditions.  Analytic  solution  shows  the  computed  field  is  accurate 
in  location  of  the  first  large  maximum,  the  amplitude  of  the  reflected  light 
and  the  amplitude  of  the  largest  maximum.  The  various  test  problems  thus 
confirm  the  validity  of  the  code. 

3.  BREAKDOWN  AND  OPTICAL  DIAGNOSTICS 

To  describe  the  nitrogen  breakdown  by  the  microwave  radiation,  we  have 


Fig.  1  —  Test  run  with  cAt  *  0.9Ay  for  *  plan*  wav*  reflected  from  a  surface  (a)  reflected 
wave  (VOS/VEO)  versus  time  (Uj)  and  (b)  electric  field  (VOS/VEO)  versus  position  (Icy) 
where  VOS  and  VEO  indicate  the  oscillatory  and  the  thermal  velocities  of  the  electron, 
respectively. 


(a)  (b) 

Fig.  2  —  Test  run  with  cAt  -  0.2Ay  for  a  plane  wave  reflected  from  a  surface  (a)  reflected 
wave  (VOS/VEO)  versus  time  (tu>)  and  (b)  electric  field  (VOS/VEO)  versus  position  (ky) 
where  VOS  and  VEO  indicate  the  oscillatory  and  the  thermal  velocities  of  the  electron, 
respectively. 


Fig.  3  —  Test  run  for  a  plane  wave  normally  incident  on  a  linear  plasma  profile. 
The  graph  shows  the  electric  field  ( VOS/VEO)  versus  position  (ky). 
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developed  an  Na  chemistry  code.  The  code  calculates  the  time  developments 
of  X*  (x),  Na+  (B) ,  N4+,  N(*D),  Na(A®E),  NaCB^r),  N»(C3tt),  Nt,  T#,  Ty,  and 

Tg,  where  the  last  three  symbols  designate  the  electron,  the  Na  -  vibrational 
and  the  gas  temperatures,  respectively.  The  purpose  for  the  calculation  of 
the  population  densities  of  several  excited  states  is  to  provide  volume 
emission  rates  for  optical  diagnostics  of  the  plasma  and  the  investigation 
of  the  ionisation  from  the  excited  states. 

3.1  Na  CHEMISTRY 

A  free  electron  gains  energy,  in  the  field  of  the  microwave  radiation, 
it  undergoes  both  elastic  and  inelastic  collisions  with  Na,  thereby  losing 
part  or  all  o£  the  energy  gained  from  the  field.  When  the  electron  energy 
reaches  and  overshoots  the  ionisation  threshold  for  Na  an  ionisation  event 
may  occur  according  to  reaction  (11) . 

e  +  Na  -4  Na+  +  2e  (II) 

Higher  energy  electrons  could  ionise  N  through  the  dissociative  ionisation. 
However,  the  threshold  for  this  process  is  higher  by  ~  10  eV  compared  to  that 
for  reaction  (11).  Therefore,  process  (11)  is  the  predominant  reaction  for 
the  breakdown  ionisation  by  microwave  radiation. 

The  deionisation  of  the  plasma  occurs  primarily  through  the  dissociative 
recombination  of  the  nitrogen  ion  Na  and  its  cluster  ion  N4+,  i.e., 

Na+  +  e  -♦  N  +  N 

N4  +  e  4  N2  +  Na 

The  N4+  cluster  forms  via  a  three-body  process 4»s  i.e., 

Na+  +  Na  +  Na  -»  N4+  +  Na  (11+) 

and  is  obviously  important  at  atmospheric  pressures  and  lower  gas  tempera¬ 
tures.  When  the  gas  temperature  rises,  however,  the  formation  of  N4+  is 
inhibited  by  the  slowing  down  of  the  forward  reaction  in  (li+)  and  by  the 
collisional  breakup  of  N4  (i.e.,  the  reverse  reaction  of  (l!+)  becomes 
operative . 


(12) 

(13) 


The  deionization  also  proceeds  via  the  three-body  collisional 
recombination,  however,  it  may  become  important  after  breakdown  and 'when 
the  electron  density  is  quite  high.  Finally,  the  electron  loss  occurs 
also  through  diffusion. 


3.2  ELECTRON  ENERGY  GAIN  AND  LOSS  PROCESSES  IN  N2 


The  electron  energy  loss  processes  in  N2  proceed  by  the  elastic  and 
inelastic  collisions.  The  energy  loss  by  elastic  collisions  is  dependent 
on  vc»  where  vc  is  the  momentum  transfer  collision  frequency.  The  collision 
frequency  depends  on  the  electron  energy  (temperature)  and  the  momentum 
transfer  cross  section.  The  cross  section6  is  shown  as  a  function  of  the 
electron  energy  in  Figure  (k) .  By  folding  this  cross  section  with  the 
appropriate  electron  velocity  distribution,  one  obtains  the  momentum  transfer 
collision  frequency.  Thus  the  rate  of  the  energy  loss  by  an  electron  through 
the  elastic  collisions  can  be  expressed  as 
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2m 
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where  m  is  the  electron  mass  and  M  is  the  mass  of  the  nitrogen  molecule. 

The  electron  inelastic  energy  loss  processes  in  N2  are  the  excitation, 
the  dissociation  and  the  ionization  of  the  nitrogen  molecule.  These 
processes  are  expressed  as 


e  +  N2 

-» 

N2  *  +  e 

(16) 

e  +  N2 

-♦ 

Na  *  +  e 

(17) 

e  +  N2 

-♦ 

N  +  N  +  e 

(13) 

in  addition  to  Equation  (1).  Equation  (16)  indicates  the  excitations  of  the 
N2  ground  state  vibrational  levels,  while  Equation  (17)  represents  the 
excitations  of  the  electronic  states  of  N2.  The  N2  electronic  states  include 
the  triplets  and  the  singlets  whose  cross  sections  are  presented  and  discussed 
elsewhere.7  The  cross  section  for  the  N2  disassociation,  however,  is  a 
composite  of  a  large  number  of  singlet  states  which  predissociate.  Thus,  the 
energy  loss  by  electrons  in  N2  can  be  summed  up  to  include  the  losses  due  to 
the  excitations  of  the  ground  state  vibrational  levels,  the  excitation  of 
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the  electronic  states  A*e,  B3e,  C3tt,  a1™,  the  dissociation,  and  the 
ionization  of  N2. 

Briefly,  the  vibrational  cross  sections  are  from  the  measurements  of 
Schulz8  and  are  well  established  as  to  the  peak  of  the  total  and  the 
shape  of  the  individual  cross  sections8  10 .  The  cross  section  for  the 
ionization  is  from  the  measurement  of  Rapp  and  Golden11  and  the  cross 
section  for  the  dissociation  is  from  the  measurements  of  Winters12  and 
Zipf  and  Mclaughlin13.  The  cross  isection  for  the  triplet  states  are 
discussed  in  the  next  section. 

The  rate  coefficients  for  all  the  inelastic  processes  discussed  above 
have  been  obtained14*  15  for  a  Maxwellian  electron  velocity  distribution  and 
are  utilized  in  the  breakdown  code.  Thus,  if  one  designates  x^  as  the 
electron  impact  rate-coefficient  for  the  excitation  of  the  ith  state  whose 
excitation  energy  is  the  rate  of  the  electron  energy  loss  in  N2  can  be 
expressed  as 


Ei  N2 


Where  the  summation  is  over  eight  ground  state  vibrational  levels,  the 
singlet  and  triplet  electronic  states,  the  dissociation  and  the  ionization 
of  N2.  A  similar  equation  for  the  energy  loss  in  the  nitrogen  atom  can  be 
utilized,  with  the  appropriate  rate  coefficients. 

Other  electron  energy  loss  processes  are  the  excitation  of  the  excited 
states  from  a  lower  state  (not  the  ground  state)  and  the  ionization  of  the 
excited  states. 

The  energy  gain  mechanisms  for  an  electron  are  the  energy  absorption 
from  the  field  of  the  microwave  i.e.,  the  inverse  bremsstrahiung,  the 
electron  impact  deexcitation  of  the  ground  state  vibrational  levels  and  the 
deexcitation  of  the  electronic  states  of  the  N2  molecule  and  the  nitrogen 
atom. 

3.3  EXCITED  STATES 

The  code  as  inuieated  earlier  calculates  the  time  histories  of  several 
excited  states  of  N2  which  are  discussed  below. 

A'^:  The  A2£  state  is  a  metastable  and  has  a  life  time  of  1  sec  and  an 
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excitation18  energy  of  6.1  eV.  Its  excitation  cross  section  by  electron 
impact  has  been  measured  and  calculated  by  several  investigators  which  are 
discussed  in  Ref  (7) »  and  shown  in  Figure  5*  The  electron  impact  ionization 
of  the  A3£  state  has  been  calculated  bu  Kukulin,  et  al.17  The  excitation 
energy  of  the  A3!  state  generally  ends  up  as  a  heating  source  for  the  gas 
because  of  the  rapid  quenching  of  the  state  by  the  atomic  nitrogen.  The 
rate  coefficient  for  the  quenching16  of  the  A3E  by  N  is  5  x  10" 11  cm3/sec. 

f B3!:  The  B3^  state  is  excited,  from  the  ground  state  of  Ng,  by  electron 
impact,  has  an  excitation  energy16  of  7.2  eV,  and  is  the  upper  level  for  the 
first  positive  band  system,  (B3rr  A3tt)  .  The  strong  vibrational  bands  of 
this  system18  are  given  in  Table  1. 

$?lC 

TABLE  1 


Relatively  Strong  Emissions  in  the  First  Positive  Band  System 


Wave  length  (A)  Transition 


7503 

6704 

6623 

6544 

646a 

6394 


Wavelength  (a)  Transition  (>6. 


4,  2 

5»  2 

6,  3 

7,  4 

8,  5 

9,  6 

10,  7 


6069 

6013 

5959 

5906 

5854 

5804 

5755 


6,  2 

7,  3 

8,  4 

9,  5 

10,  6 

11.  7 

12,  e 


The  life  times  of  these  transitions  are  in  the  range  of  10"  5  sec  to 
10"6  sec  and  are  quenched23  by  N2  with  a  rate  coefficient  of  ~  (1-2)  x  10 "11 
cm3/ sec. 

The  electron  imoact  excitation  cross  section  for  the  B3tt  state  have  been 
measured  and  calculated  by  numerous  investigators  and  some  of  these  are 
shown  in  Figure  6. 

C3t:  The  C3t  state  is  excited  from  the  ground  state  of  Ng,  by  electron 

impact  has  an  excitation  energy16  of  11  eV  and  is  the  upper  level  for  the 
second  positive  band  systems,  (C3rr  -*  33tt)  .  Some  of  the  strong18  bands  of 
this  system  are  given  in  Table  2. 
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TABLE  2 


Relatively  Strong  Emissions  in  the  Second  Positive  Band  System 
Wavelength  (*)  Transition  (■&.  Wavelength  (l)  Transition 


4059 

0,  3 

3755 

1,  3 

3998 

1,  4 

3710 

2,  4 

3943 

2,  5 

3577 

0,  1 

3894 

3,  6 

3536 

1.  2 

3804 

0,  2 

3371 

0,  0 

3159 

1,  0 

3163 

2,  1 

The  lifetimes  of  these  transitions  are29  in  the  range  of  tens  of 

nanoseconds  to  fractions  of  a  microsecond.  Those  vibrational  levels  are 

quenched  by  Na  where  the  rate  coefficient  for  the  quenching  of  the  (0,  0) 

transition  is21  1.15  x  10"11  cm3/sec. 

The  electron  impact  excitation  cross  section  for  the  second  positive 

band  system  has  been  calculated  and  measured  extensively,7  especially  the 

cross  section  for  the  excitation  of  the  (0,  0)  transition  at  3371  A.  The 

excitation  cross  section  is  shown  in  Figure  7» 

Na+(B):  The  nitrogen  molecule  has  several  ionization  continua  where 

upon  the  ionization  the  molecular  ion  is  in  an  excited  state.  These  are 
+  +  +  + 

Na  (A),  Na  (B),  Na  (C),  etc.  The  Na  (B)  state  has  an  ionization  threshold  of 

13.3  eV  and  is  the  upper  level  for  the  first  negative  band  syster 

Ns  (B2^  -*  X2!).  The  strong18  emissions  from  this  band  system  are  given  in 


Table  3. 


ERGY  (eV) 

a  of  N2(C3n ).  For  details  see  Ref. 


TABLE  3 


Relatively  Strong  Eioissions  in  the  First  Negative  Band  System 


Wave  length  (l) 

Transition  (4.  <S) 

Wavelength  (X) 

Transition  (4.  $) 

5228 

0,  3 

1*23  6 

(1,  2) 

5076 

2.  5 

3911* 

(0,  0) 

1*709 

0,  2 

1*599 

2,  1* 

1*515 

k.  6 

1*278 

0,  1 

The  lifetimes  of  these  states  range23  from  tens  of  nanoseconds  to 
microseconds  and  are  quenched  rapidly  by  Ns.  The  rate  coefficient22  for  the 
quenching  of  the  zeroth  vibrational  level  is  1*  x  10  "l0  cm3/ sec. 

The  fractional  ionization  leading  to  the  Ns  (3)  state  is  well 
established7  for  ionization  by  electrons  with  energy  of  100  eV  and  higher. 
However,  some  data  exist23  for  the  fractional  ionization  below  100  eV,  and 
emission  cross  section 24  which  should  be  of  more  relevance  to  the  ionization 
of  N2  by  the  microwave  radiation. 

N(aD):  The  metastable  state  of  the  nitrogen  atom  has  an  excitation 
energy1®  of  2.37  eV  and  arises  by  the  dissociation  of  Na,  the  dissociative 
recombination  of  Na  ,  and  by  direct  electron  impact  excitation  from  the 
ground  state  of  N.  The  electron  impact  excitation  cross  section  of  NC^) 
has  been  calculated.7 

3.1*  GAS  HEATING 

One  of  the  basic  interests  in  the  air  breakdown  is  the  heating  of  the 
plasma  and  the  neutral  species.  The  heating  of  the  neutral  species  and  the 
ions  (since  the  collision  between  heavy  particles  is  very  efficient  we 
assume  that  the  ions  and  the  neutrals  have  the  same  temperature)  occurs 
through  several  atomic  and  molecular  processes.  These  are:  the  elastic 
collisions  of  the  electrons  with  the  heavy  particles,  the  quenching  of  the 
excited  electronic  states  by  the  neutral  species,  the  dissociative  recombina¬ 
tion  of  the  molecular  ions  and  the  charge  exchange  process.  The  quenching  of 
the  vibrational  energy  through  the  collision  between  neutral  species  and 
vibrationally  excited  molecules  is  another  essential  source  for  the  heating 
of  the  gas.  However,  the  slowest  heating  rate  among  these  various  processes 
is  the  vibrational  -  translational  energy  exchange. 
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When  ch«  gas  temperature  rises,  two  additional  ionisation  processes 
become  operative.  These  are  the  thermal  ionisation  and  the  associative 
ionisation  as  expressed  by  Equations  (20)  and  (21),  respectively. 

Na  +  Na  Na+  +  Na  +  e  (20) 

N  +  N  Na+  +  e  (21) 

k.  APPLICATIONS  TO  BREAKDOWN  STUDIES 

There  exist  considerable  experimental  88 » 80  data  on  the  microwave 
breakdown  in  air  and  other  gaseous  elements.  With  the  advent  of  lasers,  the 
air  breakdown  studies  at  optical  frequencies  have  increased  considerably?7"29 
Theoretical  calculations  for  the  air,  or  gas  breakdown  generally  proceed  from 
the  Boltzman's  equation  which  is  solved  self  consistently  for  the  electron 
energy  and  its  distribution.  Analytical  methods80* 87,30  «re  also  utilised  to 
obtain  the  breakdown  threshold  power.  In  this  report  we  have  assumed  a 
maxwellian  electron  velocity  distribution  in  our  codes  (we  are  planning  to 
develop  a  Boltzman  code  for  this  purpose)  and  have  utilized  the  code  for  the 
breakdown  calculations. 

b-.l  MICROWAVE  BREAKDOWN  IN  Na 

In  order  to  compare  our  results  with  experiments,  we  have  incorporated 

the  effect  of  diffusion  into  the  code.  Thus  we  define  an  effective  ionization 

rate  «  v,  •  D  n  to  allow  for  diffusion  in  our  code.  The  diffusion  is 

assumed  to  be  free  with  a  diffusion  coefficient  D  V2  t2  =*  5  x  101S  T  . 

e  c  e 

The  time  t  *  1/v  is  the  electron-neutral  (electron- ion)  collision  time  and 

C  v 

is  the  electron  thermal  velocity.  The  diffusion  length  A  is  regarded  as 
being  specified  by  the  cavity  and  is  assumed  known.  Our  definition  of  break¬ 
down  is  defined  as  the  electric  field  required  to  achieve  a  net  increase  in 
the  electron  density  for  a  ten  nanosecond  pulse.  Once  the  cascade  process 
begins,  the  electron  density  quickly  reaches  the  critical  density.  The 
dominant  loss  processes  at  low  and  high  pressures  (Pal  Torr)  are  diffusion 
and  dissociative  recombination  (N*+  +  e  -*  Na  +  N2)  ,  respectively.  Our 
comparison  will  be  to  experimentally  measure  cw  breakdown  in  nitrogen. 

The  criterion  for  cw  breakdown  in  nitrogen  is  +  D/A8  where  is  the 

dissociative  recombination  rate31,  vf  =  5*0  x  10~3  N4+/Tg,  where  T^  is  in  eV, 
Both  the  diffusion  length  A  and  the  ambient  electron  density  (n^  *»  102  cm-3) 
were  specified  in  the  experiment  and  are  known  parameters  in  our  calculation. 


Calculations  wars  performed  for  various  frtqutncies,  diffusion  langths, 
and  prassuras.  Figura  (8)  shows  tha  rasults  obtained  by  MacDonald  at  al.*° 
for  cw  breakdown  fields  in  nitrogen  at  frequencies  992MH  and  9.4GH  with 
diffusion  lengths  of  0.631  ctr  and  0.4  cm,  respectively.  The  stars  in 
Figure  (8)  denote  the  calculated  breakdown  threshold  Intensity  using  N2-Chem. 
The  calculated  thresholds  are  higher  as  expected  for  pulsed  breakdown.  Table 
4  shows  the  electron  temperature,  breakdown  intensity,  and  characteristic 
time  for  electron  production  t  *  (dN^/dt)"1  as  a  function  of  the  back¬ 

ground  pressure  at  f  *  922  MH  and  A  *  0.631  cm. 

Table  4 

Pulse  Breakdown  Threshold  Power  for  10  cm  Wavelength  for  Various  Nitrogen 

Pressures 


Neutral^ 
densitv  (cm  3) 

Electron 
Temoerature  (eV) 

Pulse  Breakdown 
Power_  (Watts/ cm*} 

t  (sec) 

1.32  x  10iS 

1.7 

2200. 

7.1  x  10*7 

1.32  x  1017 

2.4 

80. 

4.3  x  10~7 

1.32  x  10ia 

6.4 

• 

o 

VO 

2.5  x  10"7 

1.32  x  101S 

41.0 

153. 

3.7  x  10"a 

O 


O 


The  characteristic  time  t  represents  roughly  the  maximum  growth  time  for  the 

electrons  at  the  start  of  the  avalanche  process.  Experiments  show  that  once 

the  pulse  length  tp  *  25y,s,  the  pulse  breakdown  threshold  approaches  the  cw 

breakdown  limit.  The  generation  time  is  approximately  t  *  t  log  2  and  we 

© 

require  20  -  25  generations  to  reach  electron  densities  of  order 

N@  —  101Z  cm"3.  Again  our  definition  of  the  breakdown  power  is  the  intensity 

required  to  achieve  a  net  increase  in  the  electron  density  in  10  ns. 

Therefore,  our  results  are  certainly  applicable  to  pulse  breakdown  when 

t  «  25jiS.  Since  this  is  easily  satisfied  by  our  results,  we  would  expect 
8 

the  higher  breakdown  field  as  characteristic  of  pulse  breakdown.  The  high 
electron  temperature  at  low  pressure  results  from  the  strong  diffusive  losses. 
The  agreement  at  lower  pressures  may  be  somewhat  fortuitous  since  we 
extrapolate  the  ionization  and  excitation  rates  to  energies  above  20  eV. 

Current  rate  coefficients  in  the  code  are  tabulated  for  T  uo  to  20  eV. 
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Fig.  8  —  Comparison  between  experimental  CW  N2  breakdown  threshold  power 
and  calculations  using  the  Nj  chemistry  cc  le,  CHEMN2.  Calculations  are 
indicated  by  circles  and  crosses  for  pulsed  and  CW  breakdown  fields,  respec¬ 
tively. 


However,  data  is  available  for  T#  above  this  temperature  and  will  ba 
incorporated  into  the  coda  if  desired. 
b,2  MICROWAVE  AIR  BREAKDOWN  OFF  A  REFLECTING  SURFACE 

We  have  also  coupled  the  E  and  M  code  to  a  detailed  air  chemistry 
code38  (CHMAIR)  to  study  the  breakdown  in  air.  We  have,  however,  chosen  to 
calculate  the  breakdown  off  a  reflecting  surface  to  illustrate  the  breakdown 
in  air  as  a  result  of  the  interference  between  the  incident  and  the  reflected 
waves.  In  these  calculations  we  neglect  the  hydrodynamic  effects  and  assume 
that  chemical  processes  are  more  effective  in  depleting  the  electron  density 
than  diffusion.  This  is  especially  true  in  microwave  absorption  experiments 
with  short  pulses  (t^  <  bO  ns)  and  rather  long  wavelength  radiation  conducted 
at  atmospheric  pressures.  Assuming  the  diffusion  length  A  X»  at 
atmospheric  condition  the  characteristic  electron  diffusion  time,  for  10  cm 
radiation,  is  roughly  t^  *•  (A/L^)*  tc  10”7  S8C  vhere  L^^  is  the  mean 
free  path  of  tho  electron.  Since  the  diffusion  time  tQ  »  t  “•  Mins,  our 
neglect  of  hydrodynamic  effect  would  seem  a  good  assumption. 

The  intensity  of  the  standing  wave  patterns  resulting  from  the 
interference  maximum  can  easily  be  above  the  threshold  for  air  breakdown. 

We  imagine  the  microwave  to  be  incident  on  a  reflecting  surface  with 
intensity  Ig  <  Ifi.  At  the  first  maximum  (one-quarter  wavelength  from  the 
surface)  the  field  intensity  is  I  »  b  Ig  with  air  breakdown  commencing 
rapidly  a  distance  from  the  surface  (Figure  9a) .  Figure  9b  shows  the 

rapid  increase  of  the  electron  density  as  the  avalanche  begins.  The  extreme 
sensitivity  of  the  air  breakdown  to  the  electric  field  essentially  localizes 
the  electron  density  resulting  in  sharp  scale  lengths,  L  »  (dN^/dy)** 1  <  X. 

In  this  simulation,  we  have  a  pulse  of  10cm  wavelength  radiation  incident 
normal  to  a  reflecting  surface  with  an  intensity  of  Ig  »  O.JMW/cm8  and  a 
pulse  duration  of  UOns.  The  initial  ambient  electron  density  is  taken  to  be 
n,  »  10cm"3  at  standard  temperature  and  pressure  (STP) . 

In  a  collisional  plasma,  v  »  yj,  the  cut  off  density  for  propagation  of 

c 

electromagnetic  radiation  is  not  sharply  defined  as  in  collisionless  plasmas. 
At  sufficiently  high  density  such  that  the  skin  depth  6  £  X,  the  breakdown 
region  closest  to  the  surface  is  decoupled  from  the  microwave  radiation. 

This  is  seen  in  Figure  10a,  where  we  show  the  square  of  the  electric  field  at 
the  end  of  a  -Or.s  pulse.  In  this  particular  simulation,  we  allow  for 
multiple  sources  (each  source  below  threshold)  such  that  the  net  intensity 
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Fig.  9 -(a)  Air  breakdown  off  a  reflecting  surface  at  atmospheric  conditions 
with  incident  power  I  *  0.5MW/cm2  at  X  *  10  cm.  (b)  Time  development  of 
the  electron  density  (N#/Nc). 


Fig.  10  — (a)  The  square  of  the  electric  field  at  the  end  of  the  pulse  for  t 
*  40  nsec,  (b)  Time  development  of  the  electron  density  (Ne/Nc). 
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Ig  >*  2MW/cm3  >  Ig.  The  corresponding  density  for  5  <  X  can  be  easily 
calculated  to  be  N  =■  500N  .  Figure  10b  shows  the  time  history  of  the  rise 
of  the  electron  density  to  N  ^  5OON  .  Though  air  breakdown  is  occurring 
throughout  the  volume  for  this  particular  simulation,  the  extreme  sensitivity 
of  the  air  breakdown  to  the  electric  field  again  localizes  the  region  of  the 
electron  density.  The  important  point  to  be  made  here  is  the  high  electron 
density.  As  will  be  shown  later,  the  reflection  coefficient  for  an 
idealized  sharp  electron  density  profile  is  insignificant  until  Ng  =»  lO5^  at 
X  =  10cm.  The  high  electron  density  makes  real  the  possibility  of  driving  a 

hydrodynamic  response  with  microwaves,  pending  the  inclusion  of  the 

> 

hydrodynamic  effects  and  various  pulse  shapes  and  pulse  lengths.  Early  work 
by  Wood  using  the  hydrodynamic  code  LASNEX  suggests  the  formation  of  shocks 
at  atmospheric  conditions  using  microwave  radiation. 
k.Z  REFLECTION  OF  MICROWAVE  RADIATION 


Once  the  electron  density  becomes  sufficiently  large  Ng  »  Nc>  the 
plasma  becomes  a  good  absorber  of  microwave  radiation.  However,  at  too  high 
densities,  the  plasma  becomes  a  good  reflector  of  microwave.  We  wish  to 
obtain  maximum  bounds  on  the  reflection  and  absorption  coefficient  of  the 
plasma  as  a  function  of  wavelength  and  electron  density.  We  proceed  further 
by  calculating  the  absorption  and  reflection  coefficient  for  an  idealized 
slab  density  profile  of  density  N  and  thickness  d  (See  Figure  11).  We 
consider  multiple  reflection  and  vary  the  wavelength,  density,  and  thickness 
of  the  microwave  plasma  system.  In  particular  we  examine  the  reflection  at 
two  microwave  frequency  of  3GH  and  30GH,  respectively.  The  electron-neutral 
collision  frequency  is  assumed  to  be  vc  =  5  x  1012  sec-1  which  is  typical  at 
ambient  atmospheric  conditions.  The  reflection  and  transmission  coefficients 
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Fig.  11  —  The  idealized  slab  density  profile  is  characteristic  of  systems 
with  high  electron  densities  and  steep  density  profiles. 


where  S  *  (1  -  n)/(l  +  n)  and  e  *  2ntud/c.  The  complex  refraction  index 


The  absorption  is  determined  by  the  condition  R  +  T  +  A  «  1  with  the 

thickness  measured  in  units  of  the  skin  depth  5,  6k  =  —  .  Since 

ni 

the  results  for  d  *  5  and  d  »  6  are  somewhat  similar,  we  present  the  results 
only  for  d  »  6. 

Figure  12  shows  the  reflection  and  absorption  coefficients,  as  a 
function  of  plasma  density  at  \  =  10cm  and  \  =  1cm  light,  respectively.  If 
we  assume  an  absolute  electron  density  of  «  10J'5cm“3,  we  see  the  reflection 
coefficient  is  roughly  60$  for  10cm  light  (N  =  10^  ).  Wnereas,  for  1cm 
light,  the  reflection  coefficient  is  approximately  25 $  at  Ng  =  1015cm"3  = 

lO2^.  Our  calculations  also  indicate  that  to  maintain  a  plasma  density  of 
S  =  10"5cm”2  requires  an  absorbed  power  of  ^ajjS  =**  O.SMW/craa.  Using  a  10cm 
light  would  require  an  incident  power  level  of  s  l.oMW/cm2  >  Ig. 

Obviously,  o-.ce  the  incident  power  is  above  the  breakdown  threshold,  the 


plasma  would  be  impossible  to  maintain.  However,  the  prospects  of 
maintaining  a  plasma  at  Ng  ■  1015cm“3  using  a  1cm  light  appears  feasible 
since  an  incident  intensity  of  I^nc  =■  lMW/cma  <  Ig  would  only  be  required. 
Though  the  electron  density  =*  I015cm  3  is  clearly  arbitrary,  it  is  not 
unreasonable  to  expect  the  plasma  to  become  a  good  reflector  when 

2  ^  4.2  +  2 

"p  *  vc  w  * 

5.  Conclusions 

We  have  presented  and  discussed  two  codes  developed  for  the  studies  of 
the  pulsed  microwave  interaction  with  Ng.  We  have' also  utilized  our  E  and  M 
code  coupled  to  a  chemistry  code  to  study  the  breakdown  in  air.  Our 
breakdown  threshold  power  calculations  agree  favorably  with  experimental  data 
in  Ng. 

Our  calculations  for  breakdown  near  a  surface  brings  forth  the  importance 
of  reflection  from  a  plasma  on  its  maintenance  and  heating.  However, 
additional  studies  in  this  subject  are  in  progress  to  delineate  this  point 
further,  including  the  hydrodynamic  effects. 
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APPENDIX  B 


Microwave  Energy  Coupling  in  a  Nitrogen  Breakdown  Plasma 


C.  L.  Yee,  A.  W.  Ail*,  and  W.  M.  Bollen 
Naval  Research  Laboratory,  Washington,  D.  C.  20375 

ABSTRACT 

Computer  simulations  of  microwave  coupling  to  a  nitrogen  breakdown  plasma 
have  been  performed  at  25  Torr.  Non-hydrodynamic  ionization  fronts  are 
observed  to  propagate  toward  the  radiation  source  under  a  variety  of 
circumstances.  Free  nitrogen  breakdown  simulations  in  a  spherical  system 
show  the  propagation  velocity  of  the  breakdown  wave  can  be  as  high  as  5  x 
10b  cm/sec.  An  elementary  theory  is  used  for  estimating  the  speed  of  the 
breakdown  wave  in  one  dimension.  The  results  are  in  reasonable  agreement 
with  breakdown  experiments. 


*  Plasma  Division,  Naval  Research  Laboratory. 
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1.0  Introduction 

.  The  pulsed  breakdown  in  air  and  other  gaseous  elements  has  been 
studied  extensively1  with  emphasis  on  the  threshold  power  for  breakdown 
and  its  dependence  on  the  gas  pressure  and  radiation  wavelength.  The 
hydrodynamic  effects  were  first  considered  by  Lin  and  Theofilos.2  Their 
calculation  of  the  gas  heating  showed  that  both  high  field  intensity  above 
the  breakdown  intensity  and  increasing  molecular  density  discouraged 
energy  deposition  in  air.  Though  the  preliminary  experimental  results 
confirmed  the  general  features  of  the  simple  theory,  the  observed  pressure 
waves  were  stronger  than  predicted.  Scharfman3  et  al.,  have  observed  that 
a  breakdown  plasma  prevents  the  transmission  of  microwave  energy  to  points 
beyond  the  plasma.  Experiments4  recently  performed  at  the  Naval  Research 
Laboratory  (NRL)  showed  that  the  microwaves  are  rapidly  decoupled  from  a 
target  or  plasma  surface.  This  rapid  decoupling  of  the  microwaves  would 
limit  the  amount  of  energy  that  can  be  deposited  into  the  gas.  Computer 
simulations  of  the  NRL  experiments  have  been  performed  at  a  pressure  of  25 
Torr  in  nitrogen.  Throughout  this  work,  wc  have  used  the  MINI  code 
developed  at  NRL  to  help  in  understanding  the  basic  interaction  processes 
in  the  experiments. 

2.0  The  Microwave  Nitrogen  Interaction  Code  (MINI) 

The  Mllcrowave  Nitrogen  Interaction  code  (MINI)  used  in  our 
studies  is  an  one-dimensional  multi-species,  multi-temperature,  hydro- 
dynamic,  nitrogen  chemistry  and  wave  optics  code.  The  species  followed 
dynamically  were  N2 ,  N£(X),  N£(B),  Nj,  N,  N(2D),  N2(A3E),  N2(B3tt), 

N2(C3tt),  and  the  electron  density  (ne).  The  model  calculates  the 
electron  (Te),  vibrational  (Tv)  and  gas  (Tg)  temperatures,  and  is 
space  and  time  dependent  with  microwave  absorption  and  reflection 
considered  In  the  wave  optics  mode.  The  detail  of  the  chemistry  and  the 
wave  optics  aspect  of  the  code  are  described  elsewhere5  and  will  not  be 
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repeated  here.  The  electrons  are  described  hydrodynamtcally  by  the 
equations 


Vne>  +  7,(neue)  *  Ce 


u  =  u  -  pE  -  DV(lnn)  -  DV(lnT  ) 
c  in  c  c  c  ©  c 


W  +  ?,fE.Se  +  p«“e  +  5e>  *  V  %  *  V*  *  h 


ce  1  -r  neTc  +  ~F  Wj 


where  pe  =  neTe  is  the  electron  pressure ,  3e  =  -encue  is  the 
electron  current  density,  Re  =  "menevc(^e  “  *s  the 
electron  momentum  transfer  rate,  ye  =  e/mevc  is  the  electron 
mobility,  and  De  =  Te/mevc  *s  el«ctr°n  diffusion  coefficient. 

The  electron  momentum  transfer  frequency  is  vc,  and  me  is  the  mass  of 
the  electron.  The  electron  energy  flux  is  -  -tceVTe  where  the 
electron  thermal  conductivity  is  iee  =  0ne^e^mevc  aru*  B  a 
constant.  The  terms  on  the  right  hand  side  of  the  energy  equation 
represent  energy  lost  due  to  heat  flow,  work  due  to  momentum  transfer, 
joule  heating  and  heating  due  to  elastic,  inelastic,  and  chemical 
processes.  The  term  Ce  represents  sources  and  sinks  terms  for  the 
electron  due  to  ionization,  recombination,  and  attachment.  The  Ee  term 
in  the  electron  energy  equation  includes  energy  losses  due  to  vibrational 
and  electronic  excitation  of  molecular  nitrogen.  The  Ce  and  Ee  terms 
are  discussed  more  fully  in  Ref.  (5).  The  average  velocity  of  molecular 
nitrogen,  um,  is  related  to  the  diffusion  and  mass-average-velocityfe  of 
the  heavy  particles  by  =  vm  +  u.  The  diffusion  velocity  is  vm, 
and  u  is  the  mass-average-velocity  of  the  heavy  particles. 


The  continuity  equation  for  each  of  the  heavy  particie  species 


3t(njL)  +  7» (n^)  *  Cj. 

The  momentum  equation  for  the  heavy  particles  is 


8^.(pu)  +  V » ( p uu  +  IT.)  =  p*£  - 


where  p  =  £iminii  p+  =  E^Z^en^,  and  IT.  *  ir_  +  p^I  are  the 

mass  density,  positive  charge  density,  and  pressure  tensor.  The  viscous 

portion  of  the  pressure  tensor  is  tt_,  *  Tgljnj  is  the  heavy 

particle  pressure,  and  I_  Is  the  unit  tensor.  The  heavy  particle  energy 

equation  is 

3t(eH)  +  +  +  =  ~^e*^  +  ^e*^e  +  -rolev(V  "  ev(Tg^  +  EH 

Tv 

EH  *  £lcl"lT9  *  7  Pu2 


where  is  the  total  energy  density,  Cj  is  the 

translation-rotational  specific  heat,  ev(Tv)  =  e{exp(e/Tv)  -  I}-1  is 
the  average  vibrational  energy  per  molecule,7  e  “  0.3  eV  is  the  quantum  of 
vibrational  energy  for  nitrogen,  xv  is  the  characteristic  time  for 
vibrational  relaxation,8  and  is  the  heating  terms  due  to  chemical, 
elastic,  and  inelastic  processes.  Again  a  more  detailed  discussion  of 
Eh  can  be  found  in  Ref.  (5).  The  heavy  particle  energy  flux  is  Qh  - 
-.*H^Tg,  and  «h  Is  the  thermal  conductivity  of  the  gas,  including 
translational  and  rotational  contribution.9  The  average  and  diffusion 
velocity  of  species  l  are  related  by  u^  =  vj  +  u.  The  diffusion 
velocity  of  the  ith  ion  species  is  v^  and <u  is  the  mass-average-velocity 
of  the  heavy  particles 

1  -  + 

—  ) .  m,  n,  u,  =  u  . 
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The  vibrational  energy  equation  may  be  written  approximately  as 


MVvM  *  ,,|V«IVJ.  *  *U 


where  (&  *  -«V^TV  Is  the  vibrational  energy  flux  and  *v  is  the 
vibrational  conductivity.  The  vibrational  temperature  can  be  calculated 
assuming  that  the  nitrogen  molecule  is  a  harmonic  oscillator  and  that  the 
vibrational  levels  have  a  Boltzman  distribution.  Hence,  the  vibrational 
energy  source  term  Ev  Is 
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where  Xy  is  the  excitation  rate  coefficient10  for  the  vth  vibrational 
level  obtained  from  the  experimentally  measured  cross-sections.11  The 
vibrational  energy  spacing  is  e  ■  0.3  eV  for  nitrogen. 


The  equations  solved  in  MINI  are: 

3  (n  )  +3  (n  u  )  *  C 
t  v  e'  k  K  e  e 

u  =  u  -nE-08  (inn  )  -  D  3  (inT  ) 
e  m  e  e  xv  e'  e  xv  e' 

3.(e  )+3{eu  +  pu  +  Q  )  s  R  u  +  3  E  +  E 
t  e  x  e  e  e  e  e'  e  e  e  t 

t  =  —  n  T  +  1  m  n  u2 
e  2ee  2eee 
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(2) 
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where  the  Kinematic  and  bulk  viscosity  coefficients  are  p  and  p  .  The 
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transport  coefficients  for  the  viscosity  and  thermal  conductivity  are 
taken  from  Bird.12 

The  dc  portion  of  the  electric  field  and  the  diffusion  veloci¬ 
ties  are  determined  by  demanding  the  total  current 

^T  *  3e  +  E A  =  0  (8) 

where  3$  =  ^ler? i^i  *s  current  density  of  the  ith  ion  species. 

In  steady  state,  each  of  the  ion  species  satisfies  an  equation 

Ui  =  um  -*•  -  Dj?(lnn^)  -  DjV(lnTg)  (9) 

where  the  Ion  diffusion  coefficient  and  mobility  are  related  to  the  colli¬ 
sion  frequency  vi  by  Dj  =  Tg/mjv^  and  p^  =  Zje/mjvj. 

The  mass  of  the  ith  Ion  species  is  mj.  Summing  Eq,  (9)  for  all  ion 
species  and  using  Eq.  (8)  gives  an  equation  for  the  amblpolar  electric 
field.  Once  the  amblpolar  field  is  found,  Eqs.  (9)  and  (2)  are  used  to 
solve  for  the  flow  velocities  of  the  ions  and  the  electron. 

The  Euler ian  difference  equations  of  Eqs.  (1)  -  (7)  are  solved 
using  the  algorithm  of  Rubin  and  Burstein.13  The  scheme  is  simple,  second 
order  in  accuracy,  and  is  relatively  stable  while  centering  the  dis¬ 
sipation  terms.  The  hydrodynamic  and  chemistry  portions  are  calculated 
separately  in  a  cycle,  assuming  the  flow  velocities  of  the  fluid  are 
weak.  The  hydrodynamics  Is  determined  from  the  updated  values  of  the 
chemistry.  Hence,  the  hydrodynamics  is  considered  as  a  small  correction 
to  the  chemistry. 
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Cascade  Ionization  of  Nitrogen  by  a  Microwave  Pulse 


Microwave  gas  breakdown  begins  when  a  small  number  of  priming 
electrons  acquire  sufficient  energy  from  the  electric  field  to  ionize  the 
gas.  These  priming  electrons  gain  energy  when  the  oscillatory  energy  of 
the  electrons  in  the  electric  field  is  randomized  by  collisions  into 
thermal  energy.  Once  an  electron  has  an  energy  in  excess  of  the  ioni¬ 
zation  energy  of  the  gas,  the  electron  can  easily  ionize  the  gas,  result¬ 
ing  in  the  generation  of  two  low  energy  electrons.  The  electron  avalanche 
can  be  described  by 


n  (t)  =  n  <0)exp(t/x  )  (10) 

e  e  v  B' 

where  ne(t)  is  the  electron  density  at  time  t,  and  is  the 
breakdown  time.  Physically,  td  is  the  time  required  for  an  electron 
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to  gain  the  energy  needed  to  ionize  the  gas.14  Fig.  1  shows  the  breakdown 
time,  tB(sec),  and  the  electron  temperature,  Te(eV),  for  conditions 
typical  of  microwave  breakdown  of  Nz.  A  high  radiation  field  would  pro¬ 
mote  the  electron  avalanche  by  increasing  the  energy  imparted  to  the 
electron  per  collision  with  the  molecule.  A  high  molecular  density,  how¬ 
ever,  requires  higher  intensities  for  the  avalance  breakdown.  At  moderate 
pressures  above  25  Torr,  the  dominant  energy  loss  mechanism  for  the 
electrons  is  the  excitation  of  the  electronic  and  vibrational  modes  of  the 
gas.  The  root -mean-square  electric  field  to  density  ratio,  E/nm 
(V-cm2),  does  not  uniquely  define  the  absorbed  energy  for  an  ac  electric 
field.  The  striped  region  in  Fig.  1  is  bounded  by  the  constant  intensity 
line  (I)  and  the  constant  pressure  line  (P).  The  rapid  variation  in 

T  with  increasing  intensity  allows  a  breakdown  threshold,  I  ,  to 
□  B 

be  assigned  at  ^  onset  of  breakdown.  The  ability  to  couple  microwave 
energy  to  a  breaKdown  gas  depends  sensitively  on  the  breakdown  time, 

Xg,  and  the  electron  density. 


The  efficiency  of  converting  microwave  energy  to  the  trans¬ 
lational  energy  of  the  gas  can  be  enhanced  by  coupling  the  microwaves  to  a 
preformed  plasma.  Plasma  maintenance  of  a  preformed  plasma  by  microwaves 
allows  the  microwave  deposition  to  be  localized  while  maintaining  a  high 
electron  density.  In  short  pulse  breakdown  experiments  with  pulse 
lengths,  tp  <  1  psec,  a  high  electron  density  is  a  necessary  condition 
for  rapid  heating  of  the  gas.  Plasma  maintenance  requires  the  breakdown 
time  tg  >  Tp  when  ne(0)  >  1.  If  the  local  electric  field  is  not 
carefully  matched  to  the  instantaneous  plasma  conditions,  the  gas  heating 
will  be  self-limited  by  the  formation  of  a  non-hydrodynamic  ionization 
front.  Gas  breakdown  off  a  reflecting  surface  without  a  preformed  plasma 
will  not  result  in  a  local  deposition  of  the  microwave  energy.  Since  the 
microwave  deposition  is  determined  self-conslstently  with  the  electron 
density  profile,  the  resulting  plasma  density  cannot  rapidly  heat  the  gas 
in  the  short  pulse  experiments. 

3.1  Gas  Breakdown  Near  a  Reflecting  Surface 

A  plane  wave  incident  normal  to  a  reflecting  surface  would  re¬ 
quire  less  power  for  gas  breakdown.  The  constructive  ihterference  of  the 
incident  and  reflected  microwave  results  in  the  characteristic  standing 
wave  pattern  for  the  electric  field  and  a  factor  of  four  increase  in  the 
field  intensity.  If  the  Intensity  is  above  the  breakdown  threshold, 
breakdown  of  the  gas  occurs  instantaneously  at  the  peak  nodes  of  the  elec¬ 
tric  field.  As  the  electron  density  rises,  these  "pancakes"  of  electrons 
begin  to  attenuate  the  radiation.  Pancakes  closest  to  the  reflecting  sur¬ 
face  rise  more  rapidly  early  in  the  pulse  (Fig.  2).  As  the  electron 
density  continues  to  rise,  attenuation  of  the  microwave  causes  the  pan¬ 
cakes  closest  to  the  surface  to  fail  in  density  below  that  of  pancakes  far 
from  the  surface.  Shorter  systems  with  fewer  pancakes  decouple  more 
rapidly  than  longer  systems.  However,  the  total  absorbed  energy  remains 
essentially  constant,  independent  of  the  system  size  for  a  fixed  pressure 
and  pulse  length.15 


The  finite  length  simulations  represent  the  planar  region  near 
the  surface  of  a  focused  microwave  system  where  the  distance  from  the  sur¬ 
face  is  less  than  the  radius  of  the  spotsize.  Fig.  2  shows  two  similar 
simulations  with  I  =  6.25  kW/cm2  for  an  one-and-a-half  and  three  wave¬ 
length  system  without  hydrodynamic  effects.  The  plasma  absorbs 
approximately  20%  of  the  total  microwave  energy  with  an  asymptotic 
absorption  efficiency  of  80%.  The  absorption  Is  regulated  by  the  electron 
density.  Shorter  systems  have  a  much  higher  electron  density  than  larger 
systems  and  can  more  readily  heat  the  gas  in  the  short  pulse  experiments. 
The  gas  is  not  expected  to  be  heated  more  than  the  observed  4%  above  the 
ambient  temperature  in  a  weakly  focused  system.  The  absorbed  energy  is 
distributed  over  a  larger  number  of  pancakes  with  a  very  low  electron  den¬ 
sity.  As  wiil  be  discussed  later,  the  low  electron  density  of  Fig.  2 
results  in  a  long  heating  time  for  the  gas  of  ^  6.5  psec  where 
Tg  =  Tg(dTg/dt) -1 .  The  time  required  to  heat  the  gas  to  an 
appreciable  temperature  in  a  weakly  focused  system  would  require  pulse 
lengths  of  Tp  >  6.5  psec.  However,  even  for  long  pulses,  a  significant 
heating  of  the  gas  is  never  realized.  Later  In  the  pulse,  the  pancakes 
furthest  away  from  the  reflecting  surface  continue  to  rise  in  density. 
Ultimately,  a  single  pancake  will  have  sufficient  density  to  reflect  the 
radiation.  This  pankcake  acts  as  a  secondary  reflecting  surface  and 
decouples  the  radiation  from  the  primary  metallic  surface.  Gas  breakdown 
experiments4  using  a  reflecting  surface  show  this  behavior.  Similarly, 
free  nitrogen  breakdown  simulations  without  a  surface  using  a  focused 
system  also  show  that  the  radiation  Is  decoupled  from  the  initial 
breakdown  site. 

3.2  Maintenance  of  a  Preformed  Plasma  by  Microwaves 

Microwave  coupling  to  a  pre-ionized  plasma  near  a  metallic  sur¬ 
face  show  some  general  characteristics.  Though  the  plasma  absorbs  the 
microwave  energy,  the  absorption  cannot  be  localized.  A  traveling  ioni¬ 
zation  front  is  observed  to  proceed  toward  the  radiation  source.  Fig.  3 
shows  the  time  history  of  the  ionization  front  in  the  strong 


mismatch  case  I  >  I  and  xn  >  t0.  The  critical  time  x_  is  the 
time  required  for  the  electron  avalanche  to  reach  the  critical  density  nc 
=  u)2me/4ire2  =  1013/X2(cm)  where  X  is  the  wavelength  of  the  incident 
radiation.  In  this  simulation  I  =  12.5  kW/cm2  >  Iq  =  3  kW/cm2,  Tp  =0.5 
Vi  sec,  t  £  -  0.16  ysec,  and  all  hydrodynamic  effects  are  suppresed.  The 
initial  electron  density  (pre-ionized  plasma)  was  initialized  with  a 
symmetric  profile  with  a  peak  density  of  0.1nc  and  FWHM  of  0.1X.  In  Fig. 

3,  the  peak  density  is  determined  by  Eq.  (10).  However,  the  breakdown  time, 
t0i  is  a  complicated  function  of  the  incident  and  reflected  microwaves. 

D 

It  is  clear  that  the  small  reflectivity  (R  <  10%)  is  important  in 
determining  the  propagation  speed  since  as  Fig.  3  shows,  the  separation 
distance  between  peaks  is  X/4.  The  ionization  front  proceeds  toward  the 
source  at  a  characteristic  speed  D(t)  =  X/4x^(t).  The  maximum  density 
is  decreasing  with  time  while  the  total  absorption  continues  to  increase 
with  time.  The  asymptotic  absorption  efficiency  is  approximately  95%.  The 
movement  of  the  ionization  front  limits  the  time  the  gas-plasma  systems 
samples  the  radiation  field.  The  sampling  time,  x«.»  is  determined  by 
a(cm)  =  JQS  D(t)dt  =  D(0)xs  where  D(0)  =  x/^x^O)  =  1.6  x  106  cm/sec  is 
the  early  breakdown  speed.  Estimating  the  absorption  length  as  a(cm)  *  \/2tt 
gives  a  sampling  time  of  x  *  x  .  In  the  pressure  regime  P  >  25  Torr, 

O  V 

the  gas  is  primarily  heated  through  the  quenching  of  the  electronic  states  by 
molecular  and  atomic  nitrogen.  The  fractional  change  in  the  gas. temperature 
is  Gg  =  (ATg/Tg)Max  =  xsEaeQO0nan0ea/cvnmTg  where 

Qbtg  is  the  quenching  rate  of  the  Nz  triplet  state  na  by  the  quenchant 
ng .  The  specific  heat  of  nitrogen  is  cv  and  ea  is  the  energy  level  of 
the  i\,  state.  As  an  example,  the  N2(C:it)  state  is  quenched  by  N2  with  a  rate 
coefficient16  of  1.2  x  10-11  cm3/sec.  The  fractional  change  in  the  gas 
temperature,  eg,  is  approximately  21%  from  the  ambient  temperature. 

The  region  of  the  maximum  gas  heating  is  at  the  initial  breakdown 
site  (*X/4  from  the  surface).  The  electric  field  is  less  than 
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1%  of  the  incident  time  average  electric  field  at  the  maximum  gas  tempera¬ 
ture  (Fig.  4a).  The  plasma  absorbs  approximately  70%  of  the  incident 
microwave  energy  and  has  an  electron  temperature  of  less  than  2.5  eV 
(Fig.  4b).  Maximum  gas  heating  requires  the  sampling  time  and  the  elec¬ 
tron  temperature  to  be  as  large  as  possible  without  decoupling  the 
system.  A  high  electron  temperature  populates  the  electronic  states  which 
are  rapidly  quenched  to  heat  the  gas.  Approximately  4%  and  40%  of  the  in¬ 
cident  energy  goes  into  heating  the  gas  and  exciting  the  vibrational 
states.  A  similar  simulation  with  the  peak  power  reduced  to  I  =  6.25  kW/ 
cm2  shows  no  movement  of  the  initial  ionization  region  while  increasing 
the  fractional  change  in  the  gas  temperature  to  27%.  Further  irradiance 
would  increase  the  gas  temperature.  However,  these  results  are  sensitive 
to  the  hydrodynamic  response  of  the  plasma.  Again,  a  similar  simulation 
with  hydrodynamic  effects  showed  a  reduction  in  the  gas  and  vibrational 
temperature.  The  microwave  plasma  system  again  decoupled;  however,  the 
hydrodynamic  motion  reduced  the  energy  input  into  the  gas  and  vibrational 
energy  to  0,1%  and  26%  respectively.  The  electron  density  scale  lengths 
are  sharp,  L/\  >  10~2  where  L  =  ne(dne/dx)~1 .  The  hydrodynamic  motion 
of  the  electrons  tends  to  reduce  the  electron  temperature  and  density  in 
regions  of  high  density.  The  slow  gas  flow  (<  300  cm/sec)  further  reduces 
the  electron  temperature  by  convecting  the  vibrational  energy  away  from 
regions  of  high  electron  density.  The  small  temperature  and  density 
changes  translate  to  greater  than  10%  reduction  in  the  population  of  the 
electronic  states.  Since  quenching  of  the  electronic  state  is  the  domi¬ 
nant  mechanism  for  rapidly  heating  the  gas,  the  gas  heating  is  negligibly 
small.  It  is  clear  that  the  motion  of  the  ionization  front  will  certainly 
limit  the  gas  heating. 

3.3  Gas  Breakdown  in  a  Focused  Microwave  System 

Simulations  of  a  focused  microwave  system  were  also  performed. 

If  the  field  intensity  in  a  focused  system  greatly  exceeds  the  breakdown 
threshold  at  the  focal  spot,  a  breakdown  wave  propagates  toward  the 
source.17  The  motion  of  the  microwave  plasma  absorbing  layer  is  not  a 


simple  breakdown  wave,  since  the  reflected  microwave  is  important  in 
determining  the  propagation  speed.  Fig.  5  shows  the  time  history  of  the 
electron  density  profile  during  the  course  of  a  Xp  =  0.5  ysec  pulse  in  a 
focused  system.  The  spherical  wave  equation  is  solved  for  the  electric 
field  to  model  the  convergence  of  the  field.  The  peak  incident  field  is  I 
=  10  kW/cm 2  with  a  system  aspect  ratio  of  2.5.  The  average  power  at  focus 
is  31.25  kW/cm2,  and  all  hydrodynamic  effects  are  suppressed.  The  very 
rapid  early  time  speed  of  the  absorbing  wave,  D  »  5  x  106  cm/sec,  gives  a 
sampling  time  of  x,,  *  1.7  x  10“8  sec.  However,  the  gas  continues  to 
be  heated  after  the  passage  of  the  field  by  the  de-excitation  of  the  elec¬ 
tronic  states.  The  maximum  gas  and  vibrational  temperatures  are  eg  =« 

10%  and  ev  »  260%  respectively.  The  quarter-wavelength  fine  structure 
in  the  density  profile  shows  the  effect  of  the  reflected  light.  The  time 
average  reflection  and  absorption  coefficient  are  19%  and  68%  respect¬ 
ively.  A  snapshot  of  the  electric  field  and  electron  temperature  (Fig.  6) 
shows  the  lack  of  field  penetration  (a  *  0.1X)  and  the  cold  plasma 
temperature,  Te  <  3.3  eV.  The  gas  and  vibrational  temperatures  again 
are  reduced  when  hydrodynamic  effects  are  included  in  the  simulations. 

The  gas  is  expected  to  be  heated,  at  most,  to  10%  of  the  initial  value 
during  the  course  of  a  0.5  ysec  pulse. 


The  propagation  velocity  of  the  breakdown  wave  in  a  focused 

microwave  system  can  be  calculated.  Simulations  show  the  absorption  is 

localized  at  the  head  of  the  breakdown  wave.  The  reflection  coefficient 

is  relatively  constant  after  the  initial  breakdown  (t  >  x^)»  The 

breakdown  time,  xD,  is  inversely  proportional  to  the  light  flux, 

□ 

1/x  «  (1  +  2Re£  +  £2)I,  where  the  reflection  coefficient  is  c2.  The 

□ 

Intensity  in  a  conically  focused  system  can  be  written  as 


I(X,t)  « 


(1  +  2Rec  +  c2 ) 
r2 


x  4>  ( t ) 
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where  $(t)  is  a  form  function  characterizing  the  pulse.  The  radius  of  the 
radiation  channel  at  a  position  X  on  the  channel  axis  is  r  *  rQ  +  XtanG. 
The  minuiiium  radius  of  the  waist  at  the  focai  point  is  rQ  and  0  is  the 
half-angle  of  the  radiation  column.  Following  Raizer,17  the  propagation 
velocity  of  the  breakdown  wave  is  approximately  D  =  r0/T^tan6.  Under 
the  assumptions  made,  the  velocity  D  would  be  an  upper  bound  on  the 
propagation  velocity  of  the  breakdown  wave. 

4.0  Simulation  of  the  NRL  Experiments 

We  compare  the  simulation  results  with  the  NRL  experiment.  The 
NRL  35  GHz  112  kW  focused  gyrotron  system  gives  an  average  power  at  the 
focal  spot  of  33  kW/cm2.  The  polarization  of  the  microwave  is  such  that 
the  electric  field  is  perpendicular  to  the  plasma  gradient  (S  polari¬ 
zation).  The  average  aspect  ratio  of  the  experiment  gives  an  equivalent 
simulation  aspect  ratio  of  2.5.  The  free  nitrogen  breakdown  experiments 
were  performed  at  25  Torr  with  a  100  pps  repetition  rate  at  rp  =  1.5 
ysec.  The  minimum  radius  of  the  waist  at  the  focal  spot  is  rQ  =  0.75  cm 
and  the  half-angle  of  the  radiation  column  satisfies  tanG  =  0.35.  Experi¬ 
mentally,  framing  camera  data  (Fig.  7)  gives  a  propagation  velocity  of  D  s 
4.6  x  106  cm/sec.  The  width  of  the  absorbing  wave  is  approximately  1.2 
cm.  The  measured  absolute  intensity  of  the  second  positive  band,  N2  2P 
(0-0)  at  3371  A  is  1.4  x  1021  eV/cm2-sec.  The  calculated  spatial  width  of 
the  3371  A  light  emission  is  0.8  cm  with  a  peak  intensity  of  1.8  x  1021 
eV/cm2-sec.  As  mentioned  previously,  the  calculated  speed  of  the  break¬ 
down  wave  is  5  x  106  cm/sec.  Experimentally,  the  inferred  breakdown  time 
is  tb  a  2.5  x  10-8  sec  or  =  7.4  x  10-7  sec.  This  is  to  be 
compared  to  the  calculated  breakdown  time  early  in  the  pulse,  tq  =4.5 
x  10-y  sec,  and  late  in  the  pulse,  t  -  6  x  10-8  sec.  The  simulation 

D 

results  are  felt  to  be  in  reasonable  agreement  with  the  experiment. 
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5.0  Conclusion 

*  In  conclusion,  computer  simulation  of  microwave  coupling  to  a 
weakly  ionized  nitrogen  gas  plasma  shows  excellent  absorption  of  the 
microwave  energy.  The  electronic  states  of  nitrogen  excited  by  electron 
impact  are  rapidly  quenched  by  molecular  and  atomic  nitrogen  to  heat  the 
gas  at  pressures  above  25  Torr.  However,  high  power  irradiance  can  easily 
produce  a  non-hydrodynamic  ionization  front  which  can  severely  limit  the 
heating  of  the  gas.  Microwave  gas  breakdown  experiments  with  pulse 
lengths  less  than  one  microsecond  are  expected  to  show  only  a  slight  heat¬ 
ing  of  the  gas. 
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FIGURE  CAPTIONS 


The  breakdown  time  and  the  electron  temperature  vs.  the  rms 
electric  field  to  density  ratio  E/nm. 

Gas  breakdown  near  a  reflecting  surface  showing  the  electron 
density  pancakes  for  (a)  one-and-a  half,  and  (b)  three 
wavelengths  planar  system.  The  critical  density  is  n_  =  1.5  x 
10»  cm-3  and  k  =  7.33  cm-1.  ° 

Time  history  of  the  ionization  front  in  a  planar  system  showing 
the  wave  motion  toward  the  radiation  source.  The  microwaves  are 
incident  from  the  left  onto  a  pre-ionized  plasma  near  a 
reflecting  surface  on  the  right  boundary.  The  normalization 
parameters  are  the  same  as  in  Fig.  2. 

(a)  The  time  averaged  electric  field  square  and  (b)  the  electron 
temperature  vs.  position.  The  normalization  parameters  are: 

Teg  =  0.25  eV,  k  =  7.33  cm-1 ,  and  the  arbitrary  units  of  <EZ> 

=  1  corresponds  to  an  average  intensity  of  9.4  kW/cm2. 

Time  history  of  the  ionization  front  in  a  spherical  system.  The 
microwaves  are  incident  from  the  left  boundary  with  an  intensity 
of  5  kW/cm2.  The  focused  average  power  is  31.25  kW/cm2  at  the 
right  side  boundary. 

(a)  The  time  averaged  electric  field  square  and  (b)  the  electron 
temperature  vs.  position  for  the  simulation  of  Fig.  5.  The 
normalization  parameters  are  the  same  as  in  Fig.  4. 

Framing  camera  data  showing  the  motion  of  the  ionization  front 
in  a  focused  free  nitrogen  breakdown  experiment.  The  top  and 
bottom  photographs  were  respectively  taken  0.2-  and  0.1-  ysec 
after  the  first  observable  breakdown. 
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Microwave  Energy  Coupling  to  Nitrogen  During  Breakdown 
Bollen,  C.  L.  Yee,  A.  W.  AliT,  H.  3.  Nagurney, 
and  M.  E.  Read* 


ABSTRACT 


Microwave  energy  coupling  to  nitrogen  for  breakdown  conditions 
have  been  investigated.  The  experiments  were  performed  in  25  Torr 
nitrogen  using  a  35  GHz  lys  microwave  pulse  focused  to  33  kW/cm2.  The 
plasma  created  during  breakdown  was  observed  to  travel  towards  the  source 
with  a  velocity  of  4.6  x  106  cm/s.  Spectroscopic  measurements  using  a 
nitrogen  helium  mixture  indicate  the  average  electron  density  over  a  2  cm 
cord  to  be  4.4  x  1012  cm-3,  and  the  electron  temperature  to  be  3.8  eV. 
Microwave  cutoff  was  observed  photographically  and  suggests  the  peak 
electron  density  to  be  the  collisional  critical  density.  Computer 
simulations  were  also  performed  and  their  results  are  in  accord  with  the 
experimental  results. 
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INTRODUCTION 


Microwave  breakdown  studies  of  gaseous  elements  have  been 
carried  out  extensively  over  a  wide  range  of  pressures  and  for  several 
microwave  frequencies  using  CW  and  pulsed  radiation  sources1 .  The  main 
emphasis  on  these  studies  was  the  determination  of  the  breakdown  power 
threshold,  and  its  dependence  on  the  gas  pressure  and  the  microwave 
frequency.  The  coupling  of  microwave  energy  into  the  breakdown  plasma  and 
neutral  gas  has  not  been  studied  in  detail.  The  reason  for  this  is  that 
until  recently  no  high  power  microwave  sources  have  been  available  to  per¬ 
form  such  studies.  Most  of  the  early  work  performed  on  breakdown 
thresholds  was  performed  using  high  Q  cavities  to  obtain  the  necessary 
electric  field  to  break  down  the  gas.  Once  breakdown  of  the  gas  occurred, 
the  Q  of  the  cavity  dropped  and  the  interaction  changed.  Using  the  NRL 
high  power  gyrotron  facility,  we  have  been  able  to  eliminate  the  need  for 
cavities  and  have  performed  experiments  using  a  focused  geometry2  to 
examine  the  coupling  of  microwave  energy  to  nitrogen  gas  during 
breakdown.  We  have  also  modeled  the  experiments  using  a  1-D  computer 
simulation  codea .  Simulations  were  performed  in  a  spherical  geometry 
using  a  self-consistent,  nitrogen  chemistry,  wave  optics,  microwave 
breakdown  simulation  code,  MINI. 

The  main  emphasis  of  tnis  paper  is  on  the  ionization  front 
created  during  nitrogen  breakdown  and  its  motion  and  plasma  properties  as 
observed  experimentally.  This  motion  is  not  a  movement  of  the  plasma  but 
rather  a  movement  of  where  the  strongest  ionization  of  the  gas  is  taking 
place.  We  believe  it  to  be  due  to  two  effects:  1)  delayed  ionization  and 
2)  reflection  from  the  plasma  itself.  Since  the  electron  density  is 
proportional  to  e^0,  where  I0  is  the  local  microwave  intensity, 
breakdown  occurs  faster  in  those  regions  of  high  intensity.  In  a  focused 
system  this  results  in  breakdown  beginning  at  the  focus.  The  breakdown 
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then  moves  back  into  the  lower  intensity  region  as  time  progresses.  If 
nothing  else  is  occurring,  one  would  then  observe  a  cone  of  plasma  late  in 
time.  However,  the  plasma  is  both  absorptive  and  reflective.  The 
reflection  results  in  standing  wave  patterns.  Thus  peaks  X/2  apart  are 
seen  to  grow  in  time.  The  first  peak  is  X/4  from  the  plasma  surface.  As 
the  density  increases  in  these  peaks,  they  become  reflective  and 
absorptive,  cutting  off  the  microwave  energy  to  peaks  behind  them.  Again, 
the  density  grows  as  eat*o,  and,  therefore,  those  peaks  nearest  the 
focus  initially  grow  fastest.  If  only  one  peak  were  reflecting,  we  would 
expect  X/2  structure.  However,  all  the  peaks  reflect  and,  therefore,  we 
expect  X/4,  the  distance  from  the  reflecting  surface  to  the  first  electric 
field  maximum,  to  occur  as  well.  Thus,  in  the  experiment,  we  would 
expect  a  background  plasma  due  to  delayed  ionization  with  X/2,  and  X/4 
structures  in  it  due  to  reflection.  Reflection  and  absorption  will  also 
prevent  microwave  penetration,  so  we  expect  to  see  the  plasma  disappear 
behind  the  front  of  the  ionization.  This  results  in  the  appearance  of 
motion.  This  motion  may  then  occur  as  X/2  or  X/4  jumps,  since  once  these 
structures  develop,  the  Intensity  at  these  points  is  larger,  due  to 
reflection,  than  the  background  intensity. 


This  ionization  front  motion  strongly  affects  the  coupling  of 
microwave  energy  to  the  gas.  It  determines  where  the  energy  will  be 
absorbed,  since  its  high  density  does  not  allow  microwave  penetration. 
Also,  since  it  is  moving,  it  spreads  the  heating  out  in  the  gas.  Motion 
towards  the  microwave  source  has  been  observed  by  Salisbury  and  Flint4, 
and  by  Beust  and  Ford5.  The  velocities  were  much  slower,  however,  V  »  3  x 
103  cm/s  for  Salisbury  and  Flint's  work,  and  V  *  600  cm/s  for  Beust  and 
Ford's  work.  Raizer6  advanced  a  theoretical  treatment  for  the  microwave 
discharge  propagation  in  high  pressure  air  (~1Atm)  to  explain  Beust  and 
Ford's  results  under  CW  conditions.  Other  work  by  Scharfman,  et  al7, 
observed  that  for  low  pressures  the  breakdown  plasma  became  opaque  to 
microwaves,  thus  preventing  penetration.  We  have  observed  the  formation 
of  the  front,  its  motion,  the  prevention  of  microwave  penetration,  and 
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measured  the  plasma  density  and  electron  temperature  of  the  front.  Good 
agreement  was  found  between  the  experiment  and  the  simulation  model. 

II.  EXPEhIMENT 

A.  Apparatus 

The  microwave  energy  coupling  experiments  were  performed  In  dry 
nitrogen  using  one  of  the  NRL  35  GHz,  high  power  gyrotrons8.  This  gyro- 
tron  is  capable  of  150  kW,  Ipsec  pulses  at  100  pps.  The  investigation, 
however,  was  performed  with  a  microwave  power  *  112  kW  using  S  polari¬ 
zation.  The  nitrogen  pressure  was  variable.  The  data  presented  is  at  a 
pressure  of  25  Torr  which  corresponds  to  v/m  “  0.65,  where  v  is  the 
collision  frequency  for  electron  momentum  transfer  and  o>  is  the  natural 
frequency  of  the  microwave  radiation.  The  microwave  radiation  is 
Introduced  Into  a  test  chamber  using  a  conical  horn  with  a  focusing 
dielectric  lens  In  front,  see  Figure  1.  The  lens  has  a  focal  length  of 
11.2  cm  and  a  diameter  of  7.62  cm.  The  focal  spot  is  elliptical  and  has 
an  approximate  gaussian  distribution  in  intensity  with  half  widths  of  1.5 
cm  and  1.1  cm.  The  3dB  focal  spot  area  is  approximately  1.7  cm2.  For  112 
kW  this  implies  33  kW/cm2  average  power  inside  the  3dB  spot.  Breakdown, 
without  a  reflecting  surface,  was  observed  to  cease  at  p  »  75  Torr  for  P  = 
112  kW.  Theoretically,  we  calculate  an  average  power  of  27  kW/cm2  is 
required  for  breakdown  at  75  Torr  in  nitrogen.  This  is  in  reasonable 
agreement  with  an  estimated  intensity  of  33  kW/cm2.  We  have,  however, 
ignored  the  gaussian  spread  and  assumed  an  even  distribution. 

The  breakdown  experiments  were  performed  with  and  without  a 
metal  surface  at  the  focal  point.  The  planar  metal  surface  could  be 
oriented  perpendicular,  90°,  or  with  an  angle  of  45°  to  the  incident 
radiation.  The  metal  surface  serves  to  set  up  a  standing  wave  pattern  in 
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front  of  the  surface  resulting  in  4  I0  at  the  peaks,  where  IQ  is  the 
intensity  of  the  incident  radiation.  Also,  a  thin  dielectric  could  be 
placed  off  the  metallic  surface.  This  was  used  in  determining  the  cause 
of  the  plasma  motion  observed  and  will  be  discussed  later. 

Observation  of  the  plasma  was  done  using  a  TRW  framing/streaking 
camera.  A  20  ns  exposure  time  was  used  in  the  framing  mode  and  a  2  ps  ex¬ 
posure  time  in  the  streaking  mode.  Due  to  the  amount  of  light  present, 
multiple  exposures  were  required  in  order  to  obtain  easily  visible  photo¬ 
graphs.  However,  some  single  shot  experiments  were  performed  and  showed 
that  the  results  were  reproducible  shot-to-shot.  The  electron  temperature 
and  density  were  measured  by  using  a  nitrogen-helium  gas  mixture.  In  a 
pure  nitrogen  mixture,  only  second  positive  nitrogen  bands  were 
observable.  By  introducing  helium,  the  5876A  line  could  be  used  with  the 
nitrogen  3371A  band  to  make  the  required  density  and  temperature 
measurements.  An  Oriel  model  7240  monochrometer,  f/3.7,  with  a  10A 
bandwidth  was  used.  The  monochromator  output  was  measured  using  a  1P28 
photomultiplier  with  a  base  modified  for  nanosecond  response  time. 

Absolute  Intensity  calibrations  were  performed  using  a  calibrated 
deuterium  arc  lamp  source9  for  3371A  and  a  calibrated  strip  tungsten  lamp 
source10  for  5876A. 

n 

Our  experiments  were  performed  at  25  Torr  where  we  had  approxi¬ 
mately  nine  times  the  intensity  required  to  break  down  the  nitrogen  at  the 
focal  spot.  With  the  presence  of  a  metal  surface,  a  factor  of  four  re¬ 
duction  In  the  power  required  for  breakdown  occurs  due  to  the  standing 
wave  pattern  which  is  set  up.  Due  to  the  focused  nature  of  the  experi¬ 
ment,  the  standing  wave  pattern's  intensity  drops  as  the  distance  from  the 
surface  is  increased.  For  the  45®  target,  in  fact,  there  is  only  a  small 
region  of  standing  waves  perpendicular  to  the  surface.  It  is  only  for  the 
90°  surface,  or  the  case  of  a  reflective  plasma  perpendicular  to  the 
microwave  beam,  that  a  standing  wave  pattern  back  the  length  of  the  beam 
line  exists.  A  delayed  breakdown  from  the  surface  back  towards  the  lens 


Is  then  expected,  since  the  growth  of  the  plasma  is  proportional  to 
e^o.  We  believe  that  fluctuations  in  the  reflective  plasma  surface 
and  UV  ionization  (if  any)  will  cause  a  smearing  out  of  the  standing  wave 
pattern. 


B.  Olscussion  of  the  Experimental  Results 

1.  Time  Resolved  Photographs  of  Breakdown 

The  observed  breakdown  patterns  for  the  three  cases  are  shown  in 
Figures  2  through  4.  In  Figure  2,  the  breakdown  without  a  surface  is 
shown.  A  conical  plasma  is  formed  with  its  tail  nearest  the  focal  point 
and  head  towards  the  lens.  This  indicates  that  a  large  fraction  of  the 
microwave  energy  is  concentrated  on  axis.  For  a  diffraction  limited  spot 
size  we  calculate  the  depth  of  field  (where  I  =  I0/e)  for  our  case  to  be 
3.5  cm  which  is  approximately  the  length  of  the  tail  in  the  early  break¬ 
down  region  observed  in  Figure  2.  The  total  length  of  the  breakdown 
region  is  5  cm.  We  calculate  a  length  of  3.7  cm  based  on  incident  in¬ 
tensity  decreasing  to  the  nitrogen  breakdown  threshold  as  we  move  away 
from  focus.  This  suggests  reflection  off  the  plasma  must  be  occurring. 

We  calculate  that  ■  19%  reflection  would  be  required;  the  simulation 
suggests  R  <  28%.  The  breakdown  begins  at  the  focal  spot,  but  since  the 
microwave  power  is  above  the  breakdown  threshold,  a  secondary  breakdown 
region  ahead  of  the  focus  forms  In  a  time  «  2.4  x  10“7s  after  the  initial 
breakdown.  As  the  electron  density  in  this  secondary  breakdown  region 
Increases  to  Nc  =  (<u2  +  v2  )me/(4ire2 ) ,  the  tail  (focus)  region 
disappears  due  to  the  microwaves  being  absorbed  and  reflected  by  the 
plasma  at  the  head.  Once  the  microwave  energy  is  stopped  from  reaching 
the  plasma  in  the  tail,  the  plasma  temperature  there  drops  rapidly  (x  *»  10 
ns),  and  the  emitted  light  ceases.  The  plasma  itself  also  begins  to 
decay,  but  on  a  longer  time  scale.  Late  in  time  we  see  that  a  single 
absorbing  layer  about  one  wavelength  thick  is  formed. 


Figure  3a  shows  the  framing  camera  data  for  the  perpendicular 
surface  experiment  and  Figure  3b  shows  a  streak  photograph  for  the  same 
case.  We  expect  a  standing  wave  to  be  set  up  off  the  surface  with  the 
first  maximum  X/4  from  the  surface  and  with  the  rest  of  the  maxima  spaced 
X/2  apart.  This  can  be  seen  in  Figure  3a.  Also  observed  is  a  "blooming" 
of  the  visible  plasma  near  the  target.  We  believe  this  is  an  effect  due 
to  the  focusing.  The  microwave  beam  becomes  more  Intense  near  focus. 
Therefore,  the  diameter  of  the  beam  capable  of  breaking  down  the  plasma 
increases.  Similar  to  the  no  target  case,  an  absorbing  front  travels  away 
from  the  target  towards  the  microwave  source.  From  the  displacement  in 
the  streak  photograph,  we  calculate  that  the  plasma  moves  away  from  the 
surface  with  a  velocity  a  4.6  x  106  cm/s.  It  continues  to  move  back  along 
the  microwave  beam  until  the  beam  energy  decreases  sufficiently  (due  to 
defocusing)  that  the  beating  with  any  reflected  microwaves  is  below  the 
breakdown  threshold.  This  "stagnation  point"  is  also  observed  for  the  45* 
surface.  Similar  behavior  is  also  observed  in  the  computer  simulation. 

The  45°  data  is  shown  in  Figure  4.  It  is  similar  to  the  90* 

data  except  fringes  occur  X/(2  cos  8)  apart,  8  =  45*,  with  the  first 

fringe  at  X/(4  cos  8)  from  the  surface.  Again,  an  absorbing  front  moving 
along  the  beam  axis  is  observed.  This  occurs  as  a  result  of  energy  re¬ 
flected  by  the  plasma  along  the  beam.  Since  we  are  far  above  threshold, 
we  believe  a  planar  plasma,  perpendicular  to  the  beam,  grows  near  the  tar¬ 
get  and  is  initiated  by  the  plasma  growth  at  the  standing  wave  maxima. 

Its  density  grows  (He— ^Nc)  and  it  becomes  reflective,  setting  up 
standing  waves.  Plasma  density  builds  at  these  maxima  cutting  off 
microwave  penetration,  and  the  process  repeats.  Visually,  this  appears  as 
plasma  motion.  We  would  expect  the  motion  to  Jump  in  either  X/2  or  X/4 
steps.  Some  of  the  45*  data  shows  X/2  bars  behind  the  main  moving 
plasma.  In  general,  however,  this  is  not  observed.  The  X/2  spacing  is  at 
about  the  limit  of  resolution.  Further,  we  believe  that  fluctuations  in 
the  reflective  plasma  surface  cause  a  smearing  out  of  the  standing  wave 


patter;,  and  results  In  the  more  solid  like  appearance  we  observe  for  the 
moving  plasma. 

One  other  possible  explanation  for  the  front  motion  is  that  UV 
emission  causes  ionization  to  occur  in  front  of  the  plasma.  This  new 
plasma  absorbs  microwaves,  heats  up,  and  then  the  process  repeats.  This 
would  result  in  a  uniform  smear  with  no  X/2  structure.  Our  data  is  some¬ 
what  suggestive  of  this.  To  check  this  idea  we  placed  a  thin  dielectric  a 
few  wavelengths  off  the  metallic  surface  to  stop  UV  penetration  but  allow 
microwave  transmission.  If  UV  ionization  was  the  method  of  propagation, 
we  expected  the  plasma  to  stop  moving  when  it  hit  the  dielectric.  The 
plasma  moved  through  the  dielectric  barrier.  We,  therefore,  believe  the 
motion  to  be  due  to  reflection  from  the  plasma  itself  resulting  in  ioni¬ 
zation  fronts.  However,  ionization  due  to  UV  may  be  responsible  for  some 
of  the  smearing  of  the  X/2  and  X/4  structure.  The  UV  ionization  mean  free 
path  is  «  0.05  cm  <X/2  for  our  conditions. 

For  the  45#  target  case  the  breakdown  was  examined  as  the  power 
was  decreased.  As  the  power  was  decreased,  the  ionization  front  was  ob¬ 
served  to  stagnate  closer  and  closer  to  the  target,  as  expected.  When  the 
Incoming  power  was  20  kW,  «  5.9  kW/cm2,  we  observed  only  two  fringes  and 
no  apparent  ionization  front.  At  30  kW  the  front  reappeared.  Incident 
radiation  on  a  metal  surface  for  breakdown  to  occur  at  25  Torr  is  ■  0.75 
kW/cm2 .  However,  we  are  only  twice  the  free  nitrogen  breakdown  threshold 
of  3  kW/cm?  ,  and  this  will  be  the  dominant  threshold  for  ionization  front 
travel  due  to  delayed  ionization. 

2.  Spectroscopic  Measurement  of  Density  and  Temperature 

In  order  to  measure  the  temperature  and  density  of  the 
electrons,  the  absolute  intensity  of  the  nitrogen  second  positive  band, 
3371A,  was  compared  to  the  absolute  intensity  of  the  neutral  helium  line 
at  5876A .  A  mixture  of  the  helium  and  nitrogen  was  used  to  perform  this 


experiment.  The  3371A  band  corresponds  to  the  (0,0)  transition  in  the 
second  positive  band  system  of  Nj.  We  obtained  the  electron  temperature 
*rom  the  ratio  of  the  two  lines11  and  the  electron  density  from  the 
absolute  Intensity  measurement  and  the  deduced  electron  temperature. 


For  an  optically  thin  line  and  assuming  a  steady  state  coronal 
model,  the  band  Intensity  at  337 1A  can  be  expressed  as: 
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where  Nj  is  the  density  of  the  nitrogen  molecules  in  the  ground  state, 

Xoc  is  the  emission  excitation  rate  coefficient,  Aq>0  is  the 
transition  rate  for  the  band  at  3371A  and  Aq  is  the  total  transition 
rate  for  the  excited  state,  E^q  is  the  photon  energy,  L  is  the  length 
of  the  emitting  plasma,  and  q  is  the  quenching  rate  coefficient  of  the 
excited  state  by  the  neutral  species.  A  similar  equation  can  be  written 
for  the  intensity  of  the  helium  line.  The  excitation  rate  coefficients 
for  3371A  and  5876A  were  obtained  using  the  measured  cross  sections12*13 
averaged  over  a  Maxwellian  electron  velocity  distribution.  The  transition 
rates  for  3371A  and  for  5876A  are  from  Refs.  (14)  and  (15)  respectively. 
The  quenching  rate  coefficient16  for  3371A  by  N2  is  1.15  x  10”11  cm3 /sec. 

The  measured  absolute  Intensities  of  3371A  and  5876A  in  a  mix¬ 
ture  of  15  Torr  N2  and  19  Torr  He  were  28.36  W/cm2sr  and  5.6  x  10“2 
W/crr^sr.  Using  the  above  analysis,  we  obtain  Tg  =  3.8  eV  and  Ne  =  4.4 
x  1012  cm-3  for  the  average  density  over  a  2  cm  chord.  If  quenching  of 
the  helium  line  by  nitrogen  is  included  the  electron  temperature  changes 
slightly,  Te  “  4.1  eV,  with  a  corresponding  slight  decrease  in  the 
calculated  electron  density.  This  electron  density  is  in  reasonable 
agreement  with  Langmuir  probe  data,  Ne  ”  2  x  1013  cm“3  and  simulation 
predictions,  Ne  =  3.6  x  1013  cm”3  since  these  last  two  are  peak  local 
densities.  It  Is  reasonable  to  expect  that  the  peak  density  is  at  least 
four  times  larger  then  the  average  density  since  the  density  is 


I 


proportional  to  e*ot  and  the  microwave  intensity  is  gaussian 
across  the  chord.  (The  two  centimeter  cord  is  slightly  larger  then  the 
FWHM  of  the  Intensity  gaussian.)  A  comparison  of  the  absolute  Intensity 
of  the  3371A  nitrogen  band  for  pure  N2  as  observed  experimentally  (Iexp  _ 
1.12  x  1020  eV/(cm2 -s-sr))  agreed  reasonably  well  with  the  simulation 
predictions  for  similar  conditions  (Isim  =  3.6  x  102°  eV/(s-cm2-sr)) . 

As  expected,  the  simulation  also  predicted  a  lower  electron  temperature, 

Te  »  3.25  eV,  for  the  pure  N2  case. 

III.  COMPUTER  SIMULATION 

The  breakdown  experiments  were  modeled  neglecting  hydrodynamic 
motion  and  solving  the  electric  field  in  spherical  geometry17.  Nitrogen 
breakdown  experiments  with  no  surface,  Figure  2,  show  the  microwaves 
penetrating  approximately  six  wavelengths  into  the  tenuous  plasma  early  in 
time.  The  experimental  light  emission  shows  t  e  electric  field  is 
essentially  close  to  planar  for  regions  near  the  focal  point.  However, 
for  distances  greater  than  4  cm  from  the  focal  point,  the  electric  field 
shows  an  abrupt  transition  to  a  spherical  (converging)  geometry.  The 
ratio  of  the  inner  and  outer  radius  in  a  spherical  system  defines  an 
aspect  ratio.  The  converging  geometry  of  the  experiment  defines  an  aspect 
ratio  of  3.2  in  the  converging  region  and  1.9  in  the  diffraction  limited 
focal  region  (approximately  planar).  We  modeled  this  using  an  aspect 
ratio  of  2.5  in  the  simulation.  (The  simulation  has  shown  that  for  larger 
aspect  ratios  the  behavior  is  similar,  but  the  front  does  not  travel  as 
far.)  The  simulation  used  a  0.5ps  pulse  length,  Tp,  at  an  angular 
frequency,  w ,  of  2.2  x  10!  1  sec"1  and  an  average  power  of  31.25  kW/cm2  at 
tl.e  focal  point.  The  region  modeled  was  three  wave  lengths  long  (2.57 
cm),  had  no  reflecting  surface,  and  a  nitrogen  pressure  of  25  Torr.  This 
simulation  showed  ionization  front  motion  similar  to  that  seen 
experimentally  in  Figure  2.  This  is  shown  in  Figure  5  as  a  time  history 
of  the  plasma  density  profile.  The  electric  field  drops  rapidly  beyond 
ihe  peak  electron  density  Nefnax  *  3.6  x  1013  cm'3  and  penetrates  about 


one  and  a  half  wavelengths  into  the  overdense  plasma  (Figure  6a).  The 
peak  electron  and  gas  temperatures  are  Te  *  3,5  eV  (Figure  6b)  and  Tg 
£  0.027  eV,  respectively.  Breakdown  simulations  with  hydrodynamic  effects 
show  similar  features. 

This  same  simulation  was  repeated  with  a  reflecting  surface  at 
the  focal  point.  Figure  7  shows  the  simulation  results  for  the  electron 
density.  Close  to  the  reflecting  surface,  the  radiation  is  essentially 
planar  and  exhibits  the  characteristic  standing  wave  pattern.  Electron 
density  "pancakes"  form  at  the  interference  maximum  of  the  electric 
field.  However,  at  high  power  the  electron  density  of  the  pancakes  rises 
to  Ne  Nc  at  which  time  the  radiation  is  partially  reflected  from  the 
density  pancakes.  This  shifts  the  breakdown  patterns  in  increments  of  a 
quarter-wavelength.  The  quarter-wavelength  fine  structure  of  Figure  8 
cannot  be  resolved  experimentally.  Similar  to  the  freo-nitrogen 
breakdown,  approximately  79%  of  the  incident  microwave  energy  is 
absorbed.  Early  in  time,  the  reflecting  surface  acts  to  enhance  the 
field.  However,  once  the  density  rises  to  absorb  the  radiation,  the  field 
reflects  from  the  density  pancakes.  Though  the  reflection  was  small  (R  < 
28%),  the  effects  of  the  reflection  can  also  be  seen  for  the  no  surface 
simulation  in  the  quarter-wavelength  fine  structure  of  Figure  5. 

The  ionization  front  is  observed  to  stagnate  at  about  the  same 
position  in  the  experiment  with  and  without  a  surface.  This  would  only  be 
true  if  the  target  is  essentially  decoupled  from  the  microwave  radiation 
and  the  absorption  efficiencies  are  approximately  equal.  The  simulations 
show  this  to  be  true  with  absorption  efficiencies  of  79%  and  72%  for 
breakdown  with  and  without  a  surface.  The  simulations  predict  a  front 
motion  of  V  «  5  x  106  cm/s  early  in  the  pulse  and  slowing  to  2  x  106  cm/s 
when  the  simulation  ends  In  0.5us.  This  gradual  slowing  of  the  ionization 
front  is  consistent  with  what  is  seen  experimentally  in  Figure  3b. 


IV. 


CONCLUSION 


In  conclusion,  we  have  observed  that  for  I  >  Ib>  Ionization 
fronts  form  and  move  towards  the  source  at  a  velocity  •  L6  x  106  cm/s. 
These  ionization  fronts  result  from  both  delayed  ionization  and  reflection 
of  microwave  energy  from  the  front.  The  front  density  can  become  quite 
large  and  is  seen  to  cut-off  microwave  penetration  in  about  one  wavelength 
in  late  time.  The  average  chord  plasma  density  was  spectroscopically 
measured  to  be  Ne  =  4.4  x  1012  cm-3  and  the  electron  temperature  to  be 
Te  =  3.8  eV  for  a  nitrogen-helium  mixture.  Langmuir  probe  measurements 
of  the  density  gave  reasonable  agreement  with  peak  densities  for  pure 
nitrogen  of  Ne  »  2  x  1013  cm“3 .  The  cut-off  of  microwave  penetration 
observed  photographically  suggests  an  electron  density  Ne  =  Nc  =  2.2  x 
1013  cm-3.  We  expect  the  electron  temperature  in  the  gas  mixture  to  be 
higher  since  N2  vibrational  and  low  lying  electronic  states  act  to 
decrease  Te  and  monoatomic  He  has  no  such  states.  Reasonable 
agreement  was  seen  with  the  simulation  for  the  velocity  of  the  ionization 
front,  Ne,  Te,  and  I3371.  The  simulation  predicts  small  reflection,  R 
<  28%,  and  low  gas  temperature,  Tg  <  0.027  eV. 


FIGURE  CAPTIONS 


Fig.  1  Schematic  of  nitrogen  breakdown  experiment.  Shown  as  an  insert 
is  the  3dB  focal  spot. 

Fig.  2  Framing  photographs  of  breakdown  at  25  Torr  N2  for  no  surface 
present,  I0  =  112  kW  or  33  kW/cm2  at  focus.  The  microwaves 
are  incident  from  the  left.  Frames  are  100  ns  apart  and  have  a 
20  ns  exposure  time.  Time  shown  on  right  refers  to  the  delay 
after  the  first  observable  breakdown  and  not  the  microwave 
pulse. 

Fig.  3a  Framing  photographs  at  25  Torr  N2  for  90°  metal  surface  present 
at  focus,  I0  =  112  kW.  Times  shown  on  right  refer  to  the 
delay  after  first  observable  breakdown. 

Fig.  3b  Streaking  photograph  at  25  Torr  ^  for  90°  metal  surface  present 
at  focus,  I0  =  112  kW.  Streaking  time  was  2us.  The  velocity 
of  motion  for  the  ionization  front  is  measured  to  be  4.6  x  106 
cm/s  from  the  photograph. 

Fig.  4  Framing  photographs  at  25  Torr  ^  for  45°  metal  surface,  IQ  = 

112  kW.  Times  shown  on  right  refer  to  the  delay  after  the  first 
observable  breakdown. 

Fig.  5  Time  history  of  the  ionization  front  in  a  spherical  system.  The 
microwaves  are  incident  from  the  left  with  an  intensity  of  5 
kW/cm2 .  The  focused  average  power  is  31.25  kW/cm2  at  the  right 
hand  boundary.  An  aspect  ratio  of  2.5  was  used  to  model  the 
experiment  and  Nc  is  the  collisionless  critical  density. 

Fig. 6a  Time  average  intensity  for  parameters  of  Figure  5.  Units  are  in 

a  =  (vos/Veo)2  such  that  a  =  1  corresponds  to  I  =  9.4 
kW/cnr  . 

Fig.  6b  Electron  temperature  for  parameters  of  Fig.  5  where  Te0  =  0.25 
eV. 

Fig.  7  Eiectron  density  profile  at  t  =  0.26ps  into  the  pulse.  Same 
parameters  as  in  Fig.  5  except  with  a  reflecting  surface. 
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